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<N ; Abstract 
> 

•/^ . The precision data on Z boson decays from LEP-I and SLC colliders are compared 

^V I with the predictions based on the Minimal Standard Theory. The Born approximation of 

\^ ' the theory is based on three most accurately known observables: G^ - the four fermion 

O ■ coupling constant of muon decay, mz - the mass of the Z boson, and a{mz) - the value of 

the "running fine structure constant" at the scale of mz- The electroweak loop corrections 

. ^ are expressed, in addition, in terms of the masses of higgs, rriH, of the top and bottom 

O ■ quarks, rrit and m^, and of the strong interaction constant as{mz)- The main emphasis of 

' , the review is focused on the one-electroweak-loop approximation. Two electroweak loops 

Q \ have been calculated in the literature only partly. Possible manifestations of new physics 

1-^ ' are briefly discussed. 

> 
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1 Introduction 

Z boson, electrically neutral vector boson (its spin equals 1) with mass mz — 91 GeV and 
width Tz — 2.5 GeV, Q occupies a unique place in physics. This heavy analog of the photon 
was experimentally discovered in 1983, practically simultaneously with its charged counterparts 
W^ bosons with mass mw — 80 GeV and width Tw — 2 GeV ^ . 

The discovery was crowned by Nobel Prize to Carlo Rubbia (for the bosons) and to Simon 
van der Meer (for the CERN proton-antiproton collider, which was specially constructed to 
produce W and Z bosons) 0. 

The extremely short-lived vector bosons (r = l/F ~ 10^^^ sec) were detected by their 
decays into various leptons and hadrons. The detectors, in which these decay products were 
observed, were built and operated by collaborations of physicists and engineers the largest in 
the previous history of physics. 

The discovery of W and Z bosons was a great triumph of experimental physics, but even 
more so of theoretical physics. The masses and widths of the particles, the cross-sections of their 
production turned out to be in perfect agreement with the predictions of electroweak theory of 
Sheldon Glashow, Abdus Salam and Steven Weinberg 0. The theory was so beautiful that its 
authors received the Nobel Prize already in 1979 [§], four years before its crucial confirmation. 

The electroweak theory unified two types of fundamental interactions: electromagnetic and 
weak. The theory of electromagnetic interaction - quantum electrodynamics, or QED, was cast 
in its present relativistically covariant form in the late 1940's and early 1950's and served as a 
"role model" for the relativistic field theories of two other fundamental interactions: weak and 
strong. 

The main virtue of QED was its renormalizability. Let us explain this "technical" term by 
using the example of interaction of photons with electrons. One can find systematic presentation 
in modern textbooks . In the lowest approximation of perturbation theory (the so called tree 
approximation in the language of Feynman diagrams) all electromagnetic phenomena can be 
described in terms of electric charge and mass of the electron (e, m). The small parameter of 
perturbation theory is the famous a = e'^/An ~ 1/137. 

The problem with any quantum field theory is that in higher orders of perturbation theory, 
described by Feynman graphs with loops, the integrals over momenta of virtual particles have 
ultraviolet divergences, so that all physical quantities including electric charge and mass of 
the electron themselves become infinitely large. To avoid infinities an ultraviolet cut-off A 
could be introduced. Another, more sophisticated method is to use dimensional regularization: 
to calculate the Feynman integrals in momentum space of D dimensions. These integrals 
diverge at D = 4, but are finite, proportional to 1/e in vicinity oi D = 4, where by definition 
2e = 4- D ^ (see Appendix A). 

The theory is called renormalizable, if one can get rid of this cut-off (or 1/e) by establishing 
relations between observables only. In the case of electrons and photons such basic observables 
are physical (renormalized) charge and mass of the electron. This allows to calculate higher 
order effects in a and compare the theoretical predictions with the results of high precision 
measurements of such observables as e.g. anomalous magnetic moments of electron or muon. 

The renormalizability of electrodynamics is guaranteed by the dimensionless nature of the 
coupling constant a and by conservation of electric current. 

After this short description of QED let us turn to the weak interaction. 



^Throughout the paper we use units in which h,c — 1. 



The first manifestation of the weak interaction was discovered by Henri Becquerel at the 
end of the XIX century. Later this type of radioactivity was called /3-decay. The first theory 
of /?-decay was proposed by Enrico Fermi in 1934 0. The theory was modelled after quantum 
electrodynamics with two major differences: first, instead of charge conserving, "neutral", elec- 
trical current of the type — e7Q,e + p7aP there were introduced two charge changing, "charged", 
vector currents: one for nucleons, transforming neutron into a proton, p7Q,n, another for leptons, 
transforming neutrino into electron or creating a pair: electron plus antineutrino, eja^. (Here 
e(e) denotes operator, which creates (annihilates) electron and annihilates (creates) positron. 
The symbols of other particles have similar meaning; 7q, are four Dirac matrices, a = 0, 1, 2, 3.) 

The second difference between the Fermi theory and electrodynamics was that the charged 
currents interacted locally via four fermion interaction: 

G ■ p^an ■ e7aZ/ + h.c. , (1) 

where summation over index a is implied (in this summation we use Feynman's convention: + 
for a = and — for a = 1, 2, 3); h.c. means hermitian conjugate. The coupling constant G of 
this interaction was called Fermi coupling constant. 

From simple dimensional considerations it is evident that dimension of G is (mass)~^ and 
therefore the four-fermion interaction is not renormalizable, the higher orders being divergent 
as G^A^, G^A^ ... . Why these divergent corrections still allow one to rely on the lowest 
order approximation, remained a mystery. But for many years the lowest order four-fermion 
interaction served as a successful phenomenological theory of weak interactions. 

It is in the framework of this phenomenological theory that a number of subsequent exper- 
imental discoveries were accommodated. First, it turned out that /?-decay is one of the large 
family of weak processes, involving newly discovered particles, such as pions, muons and muonic 
neutrinos, strange particles, etc. Second, it was discovered that all these processes are caused 
by selfinteraction of one weak charged current, involving leptonic and hadronic terms. Later 
on, when the quark structure of hadrons was established, the hadronic part of the current was 
expressed through corresponding quark current. Third, it was established in 1957 that all weak 
interactions violate parity conservation P and charge conjugation invariance C. This violation 
turned out to have a universal pattern: the vector form of the current V, introduced by Fermi, 
was substituted by one half of the sum of vector and axial vector. A, which meant that 7q, 
should be substituted by |7a(l + 75). 

In other words one can say that fermion ip enters the charged current only through its 
left-handed chiral component 

^L = ^(l + 75)^ (2) 

From such structure of the charged current it follows that the corresponding antifermions 
interact only through their right-handed components. 

Attempts to construct a renormalizable theory of weak interaction resulted in a unified 
theory of electromagnetic and weak interactions - the electroweak theory [0, ^ with two major 
predictions. The first prediction was the existence alongside with the charged weak current 
also of a neutral weak current. The second prediction was the existence of the vector bosons: 
Z coupled to the weak neutral current and W^ and W~ coupled to the charged current (like 
|e7Q,(l + 75)1^) and its hermitian conjugate current (like |z/7q,(1 + 75)6). 

The vector bosons were massive analogues of the photon 7; their couplings to the corre- 
sponding currents, / and g, were the analogues of the electric charge e. Thus az = f'^/4:iT 



and aw = g^/4:TT were dimensionless like a = e'^/An, which was a necessary (but not sufficient) 
condition of renormahzabihty of the weak interaction. 

The first theory, involving charged vector bosons and photon, was proposed by Oscar Klein 
just before World War II ^j. Klein based his theory on the notion of local isotopic symmetry: 
he considered isotopic doublets (p, n) and (i^, e), and isotopic triplet {B^,A^,B~). Where 5^ 
denoted what we call now W^, while A^ was electromagnetic field. He also mentioned the 
possibility to incorporate a neutral massive field C" (the analogue of our Z^). In fact this was 
the ffist attempt to construct a theory based on a non-abelian gauge symmetry, with vector 
fields playing the role of gauge fields. The gauge symmetry was essential for conservation of 
the currents. Unfortunately, Klein did not discriminate between weak and strong interaction 
and his paper was firmly forgotten. 

The non-abelian gauge theory was rediscovered in 1954 by C.N.Yang and R.Mills |^ and 
became the basis of the so-called Standard Model with its colour SU{3)c group for strong inter- 
action of quarks and gluons and SU{2)i x U{1)y group for electroweak interaction (here indices 
denote: c - colour, L - weak isospin of left-handed spinors, and Y - the weak hypercharge). 
The electric charge Q = T^ + Y/2, where T3 is the third projection of isospin. Y = 1/3 for a 
doublet of quarks, F = — 1 for a doublet of leptons. As for the right-handed spinors, they are 
isosinglets, and hence 

Y{un) = , Y{en) = -2 , Y{ur) = 4/3 , Y{dn) = -2/3 . 

Thus, parity violation and charge conjugation violation were incorporated into the foundation 
of electroweak theory. 

Out of the four fields (three of SU{2) and one of ^(1), usually denoted by W~^, W^, W" and 
B^, respectively) only two directly correspond to the observed vector bosons: W~^ and W~ . 
The Z^ boson and photon are represented by two orthogonal superpositions of W^ and B^: 

(3) 
A° = sW^ + cB^ , 

where c = cos 6, s = sin 6, while the weak angle ^ is a free parameter of electroweak theory. 
The value of 6 is determined from experimental data on Z boson coupling to neutral current. 
The "Z-charge", characterizing the coupling of Z boson to a spinor with definite helicity is 
given by 

m-Qs') , (4) 

where 

f = 9/c. (5) 

Note that the "Z-charge" is different for the right- and left-handed spinors with the same 
value of Q because they have different values of T3. The coupling constant of VT-bosons is also 
expressed in terms of e and 6: 

g = e/s . (6) 

The theory described above has many nice features, the most important of which is its 
renormahzabihty. But at first sight it looks absolutely useless: all fermions and bosons in it are 
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We denote by e,f,g the values of the corresponding charges at m^-scale, while e, /, g refer to values at 



vanishing momentum transfer. The same refers to a, az, aw ^-nd a, az, olw (see below Eqs. (mm - (18)) 



inassless. This drawback cannot be fixed by simply adding mass terms to the Lagrangian. The 
mass terms of fermions would contain both ipi and il)pi and thus explicitly break the isotopic 
invariance and hence renormalizability. The gauge invariance would also be broken by the mass 
terms of the vector bosons. All this would result in divergences of the type A^/m^, A'^/m^, etc. 
The way out of this trap is the so called Higgs mechanism [0. In the framework of the 
Minimal Standard Model the problem of mass is solved by postulating the existence of the 
doublet ifii = {ip~^,ip^) and corresponding antidoublet {i^^,—ip~) of spinless particles. These 
four bosons differ from all other particles by the form of their self-interaction, the energy 
of which is minimal when the neutral field ipi = -j^if^ + <^°) has a nonvanishing vacuum 
expectation value. The isospin of the Higgs doublet is 1/2, its hypercharge is 1. Thus, it 
interacts with all four gauge bosons. In particular, it has quartic terms \g'^WW ip\, ^pZZipl, 
which give masses to the vector bosons when ipi acquires its vacuum expectation value (vev) tj: 

mw = gri/2 , mz = fri/2 . (7) 

The magnitude of rj can easily be derived from that of the four-fermion interaction constant 
G^ in muon decay: 

G 

-^ ■ '^M7a(l + 75)/i ■ e7„(l + 75)l/e • (8) 

In the Born approximation of electroweak theory this four-fermion interaction is caused by an 
exchange of a virtual W -boson (see Figure 1). Hence ^ 

% = -^ = -^ , V = {V2G.)-'/^ = 246 GeV . (9) 

Such mechanism of appearance of masses of W and Z bosons is called spontaneous symmetry 
breaking. It preserves renormalizability |ll]]. (As a hint, one can use the symmetrical form of 



Lagrangian by not specifying the vev rj.) 

The fermion masses can be introduced also without explicitly breaking the gauge symmetry. 
In this case the mass arises from an isotopically invariant term /y ■ ^Ph'^^l'^r + h.c, where /y 
is called Yukawa coupling. The mass of a fermion m = fyrj/y^. There is a separate Yukawa 
coupling for each of the known fermions. Their largely varying values are at present free 
parameters of the theory and await for the understanding of this hierarchy. 

Let us return for a moment to the vector bosons. A massless vector boson (e.g. photon) 
has two spin degrees of freedom - two helicity states. A massive vector boson has three spin 
degrees of freedom corresponding, say, projections ±1, on its momentum. Under spontaneous 
symmetry breaking three out of four spinless states, v^^, v^^ = 771 ('^'^ ~ '^°) become third com- 
ponents of the massive vector bosons. Thus, in the Minimal Standard Model there must exist 
only one extra particle: a neutral Higgs scalar boson, or simply, higgs, Q representing a quantum 
of excitation of field (fl over its vev rj. The discovery of this particle is crucial for testing the 
correctness of MSM. 

The first successful test of electroweak theory was provided by the discovery of neutral 
currents in the interaction of neutrinos with nucleons |T^. Further study of this deep inelastic 
scattering (DIS) allowed to extract the rough value of the sin 6': s^ ~ 0.23 and thus to predict 



■'More on electroweak Born approximation, for which equality g = g holds, see below, Eqs. (14h - (pa). 
^Throughout this review we consistently use capital "H" in such terms as "Higgs mechanism" , "Higgs boson" , 
"Higgs doublet" , but low case "h" of higgs is used as a name of a particle, not a man. 



the values of mw — 80 GeV and mz — 90 GeV, which served as leading lights for the discovery 
of these particles. 

A few other neutral current interactions have been discovered and studied: neutrino-electron 
scattering |T3], parity violating electron-nucleon scattering at high energies [111], parity violation 



in atoms [15]. All of them turned out to be in agreement with electroweak theory. A predom- 



inant part of theoretical work on electroweak corrections prior to the discovery of the W and 



Z bosons was devoted to calculating the neutrino-electron [|T6[ and especially nucleon-electron 
Tn interaction cross sections. 



After the discovery of W and Z bosons it became evident that the next level in the study 
of electroweak physics must consist in precision measurements of production and decays of Z 
bosons in order to test the electroweak radiative correction. For such measurements special 
electron-positron colliders SLC (at SLAG) and LEP-I (at GERN) were constructed and started 
to operate in the fall of 1989. SLG had one intersection point of colliding beams and hence one 
detector (SLD); LEP-I had four intersection points and four detectors: ALEPH, DELPHI, L3, 
OPAL. 

In connection with the construction of LEP and SLG, a number of teams of theorists carried 
out detailed calculations of the required radiative corrections. These calculations were discussed 
and compared at special workshops and meetings. The result of this work was the publication 
of two so-called 'GERN yellow reports' []18|, [[l^, which, together with the 'yellow report' [pO| , 
became the 'must' books for experimentalists and theoreticians studying the Z boson. The 
book |2l| (which should be published in 1999) summarizes results of theoretical studies. 

More than 2.000 experimentalists and engineers and hundreds of theorists participated in 
this unique collective quest for truth! 

The sum of energies of e+ + e^ was chosen to be equal to the Z boson mass. LEP-I was 
terminated in the fall of 1995 in order to give place to LEP-II, which will operate in the same 
tunnel till 2001 with maximal energy 200 GeV. SLG continued at energy close to 91 GeV. 

The reactions, which has been studied at LEP-I and SLG may be presented in the form (see 
Figure 2) : 

e+e- ^Z^ff, (10) 

where 

// = vviveVe, v^v^, Vr^r) - invisible , 
ll{ee, fj,fl,Tf) - charged leptons , 
qq{uu, dd, ss, cc, bb) - hadrons . 

About 20,000,000 Z bosons have been detected at LEP-I and 550 thousands at SLG (but 
here electrons are polarized, which compensates the lower number of events). 

Experimental data from all five detectors were summarized and analyzed by the LEP Elec- 
troweak Working Group and the SLD Heavy Flavour and Electroweak Groups which prepared 
a special report "A Gombination of Preliminary Electroweak Measurements and Gonstraints 
on the Standard Model" [|22|. These data were analyzed in [^ by using ZFITTER code (see 
below section 6.1) and independently by J.Erler and P.Langacker [|23| . 

Fantastic precision has been reached in the measurement of the Z boson mass and width 



mz = 91,186.7(2.1) MeV, F^ = 2,493.9 ± 2.4 MeV . (11) 

Of special interest is the measurement of the width of invisible decays of Z: 

rinvisible = 500.1 ±1.9 MeV . (12) 



By comparing this number with theoretical predictions for neutrino decays it was estabhshed 
that the number of neutrinos which interact with Z boson is three [Nv = 2.994 ± 0.011). This 
is a result of fundamental importance. It means that there exist only 3 standard families (or 
generations) of leptons and quarks 0. Extra families (if they exist) must have either very heavy 
neutrinos (m^v > '^z/2), or no neutrinos at all. 

This review is devoted to the description of the theory of electroweak radiative corrections 



in Z boson decays and to their comparison with experimental data [^. Our approach to the 
theory of electroweak corrections differs somewhat from that used in [|T^-[^. We believe that 
it is simpler and more transparent (see Section 6.1). In Section 2 we introduce the basic input 
parameters of the electroweak Born approximation. In Section 3 we present phenomenological 
formulas for amplitudes, decay widths, and asymmetries of numerous decay channels of Z 
bosons. The main subject of our review is the calculation of one-electroweak loop radiative 
corrections to the Born approximation. In Section 4 they are calculated to the hadronless 
decays and the mass of the W boson, while in Section 5 - to the hadronic decays. In Section 6 
the results of the one-electroweak loop calculations are compared with the experimental data. 
Section 7 gives a sketch of two-electroweak loop corrections and of their influence on the fit of 
experimental data. Section 8 discusses possible manifestations of new physics (extra generations 
of fermions and supersymmetry). Section 9 contains conclusions. 

In order to make easier the reading of the main text, technical details and derivations are 
collected in Appendices. 

2 Basic parameters of the theory. 

The first step in the theoretical analysis is to separate genuinely electroweak effects from purely 
electromagnetic ones, such as real photons emitted by initial and final particles in reaction (|I0|) 
and virtual photons emitted and absorbed by them. The electroweak quantities extracted in 
this way are called sometimes ||2^ pseudoobservables, but for the sake of brevity we will refer 



to them as observables. 

A key role among purely electromagnetic effects is played by a phenomenon which is called 
the running of electromagnetic coupling "constant" a(5'^). The dependence of the electric 
charge on the square of the four-momentum transfer q^ is caused by the photon polarization of 
vacuum, i.e. by loops of charged leptons and quarks (hadrons) (see Figure 3(a)). 

As is well known (see e.g. 



a = a{q^ = 0) = [137.035985(61)]"^ . (13) 

It has a very high accuracy and is very important in the theory of electromagnetic processes 
at low energies. As for electroweak processes in general and Z decays in particular, they are 
determined by p2 

a = a(g2 = m\) = [128.878(90)]"^ , (14) 

the accuracy of which is much worse. 
It is convenient to denote 

^Combining eq.(O) with the data on i^^e^ p3and Ve^" p^ scattering allowed to establish that z/g, t'p and 
Vt have equal values of couplings with Z boson |2q] . 



where 

5a = 6ai + 5ah = 0.031498(0) + 0.05940(66) (16) 



(for the value of 6ai see ref . ||28| , while for the value of 6ah see p9[ ) . 

It is obvious that the uncertainty of 6a and hence of a stems from that of hadronic contri- 
bution 6ah- 

While a(g^) is running electromagnetically fast, aziq"^) and awiq"^) are "crawling" elec- 
troweakly slow for g^ < ?7i|: 

(17) 

(18) 

The small differences az — «z and aw — o^w are caused by electroweak radiative corrections. 
Therefore one could and should neglect them when defining the electroweak Born approxima- 
tion. (We used this recipe {g = g) when deriving the relation (||) between G^ and rj.) 

The theoretical analysis of electroweak effects in this article is based on three most accurately 
known parameters G^, a (Eq.(|l^)) and mz (Eq.(|n|)). 

G^ = 1.16639(1) ■ 10"^ GeV-2 . (19) 



az = «z(0) -- 


= 1/23.10 


az = aziml) = 


= 1/22.91 


aw = aw{0) = 


: 1/29.01 , 


aw = aw{ml) = 


= 1/28.74 



This value of Gn [^ is extracted from the muon life-time after taking into account the 



purely electromagnetic corrections (including bremsstrahlung) and kinematical factors |^ 

25 



h^"^ 

> 


192773 ^^m2/^ 277 ^ 


m 


= 1 - 8x + 8x^ - x^ - Ux^logx 


a{m^)~ 


1 -1 2 TTlf, 1 

= a log '^ + ^ 1 

377 rUp 677 



4 )] , (20) 



where 
and 



Now we are ready to express the weak angle 6 in terms of G^, a and mz- Starting from 
Eqs.(|), (^ and (|^), we get in the electroweak Born approximation: 



G = = = ('21) 

" 4v/2m2 4v^m2 V2m%s^c^ ^ ' 



9^ _ P 77a 

Wherefrom: 

f = AV2G^ml = 0.54863(3) , (22) 

7=0.74070(2) 

sin^ 29 = 4:na/V2G^ml = 0.71090(50) , (23) 

s^ = 0.23116(23) , (24) 

c = 0.87683(13) . (25) 



The angle 9 was introduced in mid-1980s |^. However its consistent use began only in the 
1990s 1^1 . Using 9 automatically takes into account the running of Q;(q'^) and makes it possible 



to concentrate on genuinely electroweak corrections as will be demonstrated below. 



( In this review we consistently use mz as defined by EWWG in accord with our eq. ( p.7| ). 
Note that a different definition of the Z boson mass rriz is known in the hterature, related to 
a different parametrization of the shape of the Z boson peak [Q . ) 

The introduction of the Born approximation described above differs from the traditional 
approach in which a — a is treated as the largest electroweak correction, masses my/ and mz 
are handled on an equal footing, and the angle 6w, defined by 

cw = cosdw = mw/mz, s^ = 1 - c^ , (26) 

is considered as one of the basic parameters of the theory. (Note that the experimental accuracy 
of 9w is much worse than that of 6*.) 

After discussing our approach and its main parameters we are prepared to consider various 
decays of Z bosons. 

3 Amplitudes, widths, and asymmetries. 

Phenomenologically the amplitude of the Z boson decay into a fermion-antifermion pair // 
can be presented in the form: 

M{Z ^ //) = i/^^(^v'/7a + ^A/7a75)^/^« , (27) 

where coefficient / is given by Eq.(p2D.^ In the case of neutrino decay channel there is no final 
state interaction or bremsstrahlung. Therefore the width into any pair of neutrinos is given by 

r, = T{Z ^ vv) = ATo{g\, + gl,) = 9>T,gl , (28) 

where neutrino masses are assumed to be negligible, and Fq is the so-called standard width: 

To = !^ = -^^ = 82.940(6) MeV . (29) 

For decays to any of the pairs of charged leptons II we have: 



Yi = T{Z -^ II) = 4Fo 



2 



n , 3a, r, , 3a; m, , 
gli{l + — + g\i{l + -r- 6^- 



(30) 



'■z 

The QED "radiator" (l + 3a/47r) is due to bremsstrahlung of real photons and emission and 
absorption of virtual photons by / and I. Note that it is expressed not through a, but through 
a. 

For the decays to any of the five pairs of quarks qq we have 

F, = T{Z -. qq) = UToiglgR^g + gl^Ry,] • (31) 

Here an extra factor of 3 in comparison with leptons takes into account the three colours of 
each quark. The radiators Raq and Rvq contain contributions from the final state gluons and 
photons. In the crudest approximation 

a 
Rvq = Raq = 1 H , (32) 



^ Z boson couplings are diagonal in flavor unlike that of W boson, where Cabibbo-Kobayashi-Maskawa mixing 
matrix pq | should be accounted for in the case of couplings with quarks. 
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where asil"^) is the QCD running couphng constant: 

as = as{q^ = m|) ~ 0.12 . (33) 

(For additional details on d^ and radiators see Appendix E.) 

The full hadron width (to the accuracy of very small corrections, see Sections 5 and 7) is 
the sum of widths of five quark channels: 

T;, = r„ + r, + r, + r, + r, . (34) 

The total width of the Z boson: 

Tz = Th + T, + T^ + T^ + 3T, . (35) 

The cross section of annihilation of e^e^ into hadrons at the Z peak is given by the Breit- 
Wigner formula 

i2vrr,r. 

Finally the following notations for the ratio of partial widths are widely used: 

Rb = — , Rc = — , Ri = — . (37) 

i- h i- h i « 

(Note that F; in the numerator of Ri refers to a single charged lepton channel, whose lepton 
mass is neglected.) 

Parity violating interference of g^f and gvf leads to a number of effects: forward-backward 
asymmetries ApB, longitudinal polarization of r- lepton Pr, dependence of the total cross-section 
at Z-peak on the longitudinal polarization of the initial electron beam Alr, etc. Let us define 
for the channels of charged lepton and light quark {u, d, s, c) whose mass may be neglected the 
quantity 

= l^Afgyf (3g) 

9Af + 9vf 
For f = b: 

A = '^9Ab9vb , . 

where f ^ is the velocity of the 6-quark: 



Vb 



\ 



1-^. (40) 

mi 



Here rhb is the value of the running mass of the 6-quark at scale raz calculated in MS 



scheme [36 



The forward-backward charge asymmetry in the decay to // equals: 

^FB - ^^ + ^^ - l^eAf , (41) 

where Nf{Nb) - the number of events with / going into forward (backward) hemisphere; A^. 
refers to the creation of Z boson in e"'"e~-annihilation, while Af refers to its decay in //. 
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The longitudinal polarization of the r-lepton in the decay Z ^ rf is Pr = —A^. If, however, 
the polarization is measured as a function of the angle 9 between the momentum of a r~ and 
the direction of the electron beam, this allows the determination of not only Ar-, but Ae as well: 

AJl + cos^d) + 2AeCose 



Pr [COS ( 



(42) 



1 + C0S2 9 + 2ArAe 

The polarization Pr is found from Pt-{cos6) by separately integrating the numerator and 
the denominator in Eq.(^) over the total solid angle. 

The relative difference between total cross-section at the Z-peak for the left- and right- 
polarized electrons that collide with non-polarized positrons (measured at the SLC collider) 
is 

A,n - ^^^ = ^ . (43) 

The measurement of parity violating effects allows one to determine experimentally the ratios 
gvf/gAf, while the measurements of leptonic and hadronic widths allow to find gAf and d^. 



Table 1 



Observable 


Experiment 


Born 


Pull 


mw [GeV] 


80.390(64) 


79.956(12) 


6.8 


mw/mz 


0.8816(7) 


0.8768(1) 


6.8 




0.2228(12) 


0.2312(2) 


-6.8 


Ti [MeV] 


83.90(10) 


83.57(1) 


3.3 


9ai 


-0.5010(3) 


-0.5000(0) 


-3.3 


9vi/9ai 


0.0749(9) 


0.0754(9) 


-0.5 


^2 


0.2313(2) 


0.2312(2) 


0.5 



Table 1 compares the experimental and the Born values of the so-called "hadronless" ob- 
servables rriyy, g^i and gyi- For the reader's convenience the table lists different representations 
of the same observable known in the literature 



o2 



1- 



ml 



'eff 



sm 



2 nlept 
^eff 



1 



1- 



9ai' 



(44) 



(45) 



The experimental values in the table are taken from ref. [Q, assuming that lepton univer- 
sality holds. The pull shown in the last column is obtained by dividing the difference "Exp 
- Born" by experimental uncertainty (shown in brackets). One can see that the discrepancy 
between experimental data and Born values are very large for mw and substantial for g^i- That 
means that electroweak radiative corrections are essential. As for gvi/gAi, its experimental and 
Born values coincide. Moreover the theoretical uncertainty is the same as the experimental one; 
thus the pull is vanishing practically. Such high experimental accuracy for gvi/gAi has been 
achieved only recently. As for m,^ and F;, their experimental uncertainties are much larger 
than the theoretical ones. 

We would like to mention that in 1991, when we published our first paper on electroweak 
corrections to Z-decays, the LEP experimental data were in perfect agreement with Born 
predictions of Tablel. This demonstrates the remarkable progress in experimental accuracy. 
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4 One-loop corrections to hadronless observables. 

4.1 Four types of Feynman diagrams. 

Four types of Feynman diagrams contribute to electroweak corrections for the observables of 
interest to us here, raw /f^z, dAi, Qvi/ dAi'- 

1. Self-energy loops for W and Z bosons with virtual v,l,q,H,W and Z in loops (see Figures 
3(b)-3(n)). 

2. Loops of charged particles that result in transition of Z boson into a virtual photon (see 
Figures 3(o)-3(r)). 

3. Vertex triangles with virtual leptons and a virtual W or Z boson (see Figures 4(a)-4(c)). 

4. Electroweak corrections to lepton and Z boson wavefunctions (see Figures 4(d)-4(g)). 

It must be emphasized that Z boson self-energy loops contribute not only to the mass mz and, 
consequently, to the mw/'^z ratio but also to the Z boson decay to //, to which Z ^^ 7 tran- 
sitions also contribute because these diagrams give corrections to the Z-boson wavefunction. 
Moreover, there is no simple one-to-one correspondence between Feynman diagrams and am- 
plitudes. This is caused by the choice of G^ as an input observable which enters the expression 
for s and c. As a result e.g. there is a contribution to niw/rnz coming from the box and vertex 
diagrams in the one-loop amplitude of the muon decay. In a similar way the self-energy of the 
W boson enters the amplitudes for decay Z -^ II. More on that see in Appendix D. 

Obviously, electroweak corrections to mw/mz,gAi and QviIqai are dimensionless and thus 
can be expressed in terms of a, c, s and the dimensionless parameters 

t={^f , h = C^f , (46) 

mz mz 

where rrit is the mass of the t-quark and itlh is the higgs mass. (Masses of leptons and all 
quarks except t give only very small corrections.) 

4.2 The asymptotic limit at rnl ^ '^^• 

Starting from papers by Veltman [^ it became clear that in the limit t ^ 1 electroweak 



radiative corrections are dominated by terms proportional to t. These terms stem from the 
violation of weak isotopic invariance by large difference of mt and mb (see Figures 3(c), 3(i), 
and 3(j)). 

After the discovery of the top quark it turned out that experimentally t ~ 3.7. As we will 
demonstrate in this review, for such value of t the contributions of the terms, which are not 
enhanced by the factor t are comparable to the enhanced ones. Still it is convenient to split the 
calculation of corrections into a number of stages and begin with calculating the asymptotic 
limit for t ^ 1. 

The main contribution comes from diagrams that contain t- and 6-quarks because large 
difference of rrit and mi, strongly breaks isotopic invariance. A simple calculation (see Appendix 
D) gives the following result for the sum of the Born and one-loop terms: 

-W-z = c+^^^^^,^J-^at . (47) 
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3 

R = 9vi/9Ai = 1 - 4s2 + —— -at , (49) 

47r(c"^ — s^j 

1 3 _ 

2 o47rs^c"' 

The presence of t-enhanced terms in radiative corrections to Z boson decays allowed to 
predict the mass of the top quark before its actual discovery [^ ^. 



4.3 The functions Ki(t, /i), VA{t,h) and VR{t,h). 

If we now switch from the asymptotic case of t ^ 1 to the realistic value of t, then one should 
make the substitution in the eqs.(|4^) - (|50|): 

t^t + Ti{t) , (51) 

in which the index i = m,A,R, v denotes mw /niz, Qau R = 9vi/ 9ai and g^, respectively. 

The functions Tj are relatively simple combinations of algebraic and logarithmic functions. 
Their numerical values for a range of values of mt are given in Table 2. The functions Ti{t) 
thus describe the contribution of the quark doublet t, 6 to my/lmz, Qa, R = QviIqai and g^- 
If, however, we now take into account the contributions of the remaining virtual particles, then 
the result can be given in the form 

t^Vi{t,h) = t + %(t) + Hi{h) + Ci + 5Vi(t) . (52) 

Here Hi{h) contain the contribution of the virtual vector and higgs bosons W, Z and H and 
are functions of the higgs mass mn- (The mass of the H^-boson enters Hi{h) via the parameter 



c, defined by equation (|25D). The explicit form of the functions H^ is given in refs. PTl H3] and 



their numerical values for various values of m^ are given in Table 3. 
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Table 2 



rat 


t 


J- m 


Ta 


Tr 


(GeV) 










120 


1.732 


0.323 


0.465 


0.111 


130 


2.032 


0.418 


0.470 


0.154 


140 


2.357 


0.503 


0.473 


0.193 


150 


2.706 


0.579 


0.476 


0.228 


160 


3.079 


0.649 


0.478 


0.261 


170 


3.476 


0.713 


0.480 


0.291 


180 


3.896 


0.772 


0.481 


0.319 


190 


4.341 


0.828 


0.483 


0.345 


200 


4.810 


0.880 


0.484 


0.370 


210 


5.303 


0.929 


0.485 


0.393 


220 


5.821 


0.975 


0.485 


0.415 


230 


6.362 


1.019 


0.486 


0.436 


240 


6.927 


1.061 


0.487 


0.456 



Table 3 



rriH 


h 


Hm 


Ha 


Hr 


(GeV) 










0.01 


0.000 


1.120 


-8.716 


1.359 


0.10 


0.000 


1.119 


-5.654 


1.354 


1.00 


0.000 


1.103 


-2.652 


1.315 


10.00 


0.012 


0.980 


-0.133 


1.016 


50.00 


0.301 


0.661 


0.645 


0.360 


100.00 


1.203 


0.433 


0.653 


-0.022 


150.00 


2.706 


0.275 


0.588 


-0.258 


200.00 


4.810 


0.151 


0.518 


-0.430 


250.00 


7.516 


0.050 


0.452 


-0.566 


300.00 


10.823 


-0.037 


0.392 


-0.679 


350.00 


14.732 


-0.112 


0.338 


-0.776 


400.00 


19.241 


-0.178 


0.289 


-0.860 


450.00 


24.352 


-0.238 


0.244 


-0.936 


500.00 


30.065 


-0.292 


0.202 


-1.004 


550.00 


36.378 


-0.341 


0.164 


-1.065 


600.00 


43.293 


-0.387 


0.128 


-1.122 


650.00 


50.809 


-0.429 


0.095 


-1.175 


700.00 


58.927 


-0.469 


0.064 


-1.223 


750.00 


67.646 


-0.506 


0.035 


-1.269 


800.00 


76.966 


-0.540 


0.007 


-1.311 


850.00 


86.887 


-0.573 


-0.019 


-1.352 


900.00 


97.410 


-0.604 


-0.044 


-1.390 


950.00 


108.534 


-0.633 


-0.067 


-1.426 


1000.00 


120.259 


-0.661 


-0.090 


-1.460 
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L'm — 


-1.3497 + 4.13(55^ 


Ca = 


-2.2621 - 2.63(55^ 


Cr = 


-3.5045 - 5.725s'^ 


a = 


-1.1641 -4.88^5^ 



The constants Ci in eq. ( [5^ ) include the contributions of hght fermions to the self-energy of the 
W and Z bosons, and also to the Feynman diagrams, describing the electroweak corrections to 
the muon decay, as well as triangle diagrams, describing the Z boson decay. The constants Ci 
are relatively complicated functions of s^ (see PT| , 0). We list here their numerical values for 
s2 = 0.23110- 6 s^: 

(53) 

(54) 

(55) 

(56) 

4.4 Corrections dVi{t). 

Finally, the last term in equation (|52D includes the sum of corrections of three different types. 
Their common feature is that they do not contain more than one electroweak loop. 

6V, = 6iV, + 62V, + 63V, (57) 

1. The corrections SiVi are extremely small. They contain contributions of the W^-boson 
and the t-quark to the polarization of the electromagnetic vacuum Swc^ and 6ta, which 
traditionally are not included into the running of <y{q'^), i.e. into a (see Figure 5). It 
is reasonable to treat them as electroweak corrections. This is especially true for the 
VT-loop that depends on the gauge chosen for the description of the W and Z bosons. 
Only after this loop is taken into account, the resultant electroweak corrections become 
gauge- invariant, as it should indeed be for physical observables. Here and hereafter in 
the calculations the 't Hooft-Feynman gauge is used. 

16 1 

SiVmit, h) = ns^-{Swa + Sta) = -0.055 , (58) 

3 a 

16 1 

6iVR{t, h) = 7Ts^c^-{6wa + 6ta) = -0.181 , (59) 

3 a 

6iVA{t,h)=6iV,{t,h) = 0, (60) 

where 

5wa = 0.00050 , (61) 

Sta - -0.00005(1) . (62) 

(See eqs.(B.18) and (B.17) from Appendix B. Unless specified otherwise, we use nit = 175 
GeV in numerical evaluations.) 

2. The corrections 62Vi are the largest ones. They are caused in the order a&s by virtual 
gluons in electroweak loops of light quarks q = u,d, s, c, b and heavy quark t (see Figure 
6): 

62V,{t)=6lV, + 6lV,{t) . (63) 



16 



Due to asymptotic freedom of QCD [^ these corrections were calculated in perturbation 
theory. The analytical expressions for corrections ^f^ and SlVi{t) are given in |^, ^. 
Here we only give numerical estimates for them, 



<5?K. 



-0.377- 



a. 



TX 



SWa = 1.750 



a.c 



TT 



5IVr = 






-11.67 



asimt 



40.10 



-11.88 



7r 


^as{mt) 


71 


as{mt) 



-10.61- 



a.« 



TT 



-9.18- 



.a.« 



-10.80- 



TT 



vr 

TT 



where EC 



a« 



aJmtj 



1 + jf^d.logt 



(64) 

(65) 
(66) 
(67) 

(68) 

(69) 

(70) 



(For numerical evaluation, we use dg = as{mz) = 0.120.) We have mentioned already 
that the corrections 5\Vi{t), whose numerical values were given in (|67D-(|69D, are much 
larger than all other terms included in 5Vi. We emphasize that the term in 5\Vi that is 



leading for high t is universal: it is independent of i. As shown in [^, this leading term 
is obtained by multiplying the Veltman asymptotics t by a factor 



27r^-f 6ds(mi) 



or, numerically. 
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t^ til -2m 



TT 



asimt) 



TC 



(71) 



(72) 



Qualitatively the factor (|7ll) corresponds to the fact that the running mass of the t- 
quark at momenta p^ ~ m^ that circulate in the t-quark loop is lower than the "on the 
mass-shell" mass of the t-quark. It is interesting to compare the correction (|7^) with the 
quantity 



m^ 



m^(p; 



-m. 



m^(l-2.7i 



,as{mt). 



TC 



(73) 



calculated in the Landau gauge in |Q, p. 102. The agreement is overwhelming. There 
is, therefore, a simple mnemonic rule for evaluating the main gluon corrections for the 
t-loop. 

3. Corrections S^Vi of the order of aa"^ are extremely small. They were calculated in the 
literature ESI for the term leading in t (i.e. aa'^t). They are independent of i (in numerical 
estimates we use for the number of light quark flavors Nf = 5): 



SsViit) ~ -(2.38 - 0.l8Nf)al{mt)t ~ -lA8al{mt)t = -0.07 



(74) 
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4.5 Accidental (?) compensation and the mass of the t-quark. 

Now that we have expressions for all terms in eq. (^) , it will be convenient to analyze their roles 
and the general behaviour of the functions Vi{t,h). As functions of rrit at three fixed values 
of iTiH, they are shown in Figures 7, 8, 9. On all these figures, we see a cusp at rrit = rriz/^. 
This is a typical threshold singularity that arises when the channel Z — i> tt is opened. It is of 
no practical significance since experiments give mt — 175 GeV. What really impresses is that 
the function Vr vanishes at this value of rrtt. This happens because of the compensation of the 
leading term t and the rest of the terms which produce a negative aggregate contribution, the 
main negative contribution coming from the light fermions (see eq. (p5D for the constant Cr). 

In the one-electroweak loop approximation each function Vi(t, h) is a sum of two functions 
one of which is t-dependent but independent of /i, while the other is /i-dependent but inde- 
pendent of t (plus, of course, a constant which is independent of both t and h). Therefore the 
curves for mn = 100 and 800 GeV in Figures 7, 8, 9 are produced by the parallel transfer of 
the curve for mn = 200 GeV. 

We see in Figure 9 that if the t-quark were light, radiative corrections would be large and 
negative, and if it were very heavy, they would be large and positive. This looks like a conspiracy 
of the observable mass of the t-quark and all other parameters of the electroweak theory, as a 
result of which the electroweak correction Vr becomes anomalously small. 

One should specially note the dashed parabola in Figures 7, 8 and 9 corresponding to the 
Veltman term t. We see that in the interval < m.^ < 250 GeV it lies much higher than Va 
and Vr and approaches Vm only in the right-hand side of Figure 7. Therefore, the so-called 
non-leading 'small' corrections that were typically replaced with ellipses in standard texts, are 
found to be comparable with the leading term t. 

A glance at Figure 9 readily explains how the experimental analysis of electroweak cor- 
rections allowed, despite their smallness, a prediction, within the framework of the minimal 
standard model, of the t-quark mass. Even when the experimental accuracy of LEP-I and SLC 
experiments was not sufficient for detecting electroweak corrections, it was sufficient for estab- 
lishing the t-quark mass using the points at which the curves VR{mt) intersect the horizontal 
line corresponding to the experimental value of Vr and the parallel to it thin lines that show 
the band of one standard deviation. The accuracy in determining rrit is imposed by the band 
width and the slope oi VR{mt). 

The dependence Vi^mn) for three fixed values of rrit = 150, 175 and 200 GeV can be 
presented in a similar manner. As follows from the explicit form of the terms Hi{mH), the 
dependence Vi{mH) is considerably less steep (it is logarithmic). This is the reason why the 
prediction of the higgs mass extracted from electroweak corrections has such a high uncertainty. 
The accuracy of prediction of mn greatly depends on the value of the t-quark mass. If rrit = 
150 ± 5 GeV, then mn < 80 GeV at the 3a level. If rrit = 200 ± 5 GeV, then rnn > 150 GeV 
at the 3(7 level. If, however, rrit = 175 ± 5 GeV, as given by FNAL experiments [^, we are 
hugely unlucky: the constraint on mn is rather mild (see Figure 10). 

Before starting a discussion of hadronic decays of the Z boson, let us 'go back to the roots' 
and recall how the equations for Vi{mt,mH) were derived. 
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4.6 How to calculate V^? 'Five steps'. 

An attentive reader should have aheady come up with the question: what makes the amphtudes 
of the lepton decays of the Z boson in the one-loop approximation depend on the self-energy 
of the W boson? Indeed, the loops describing the self-energy of the W boson appear in the 
decay diagrams of the Z boson only beginning with the two-loop approximation. The answer 
to this question was already given at the beginning of Section 4. We have already emphasized 
that we find expressions for radiative corrections to Z-boson decays in terms of a, mz and G^. 
However, the expression for G^ includes the self-energy of the W boson even in the one-loop 
approximation. The point is that we express some observables (in this particular case, mw/fnzi 
9 Ah 9vi/ 9ai) in terms of other, more accurately measured observables (a, m^, G^). 

Let us trace how this is achieved, step by step. There are altogether 'five steps to happiness', 
based on the one-loop approximation. All necessary formulas can be found in Appendix D. 

Step I. We begin with the electroweak Lagrangian after it had undergone the spontaneous 
violation of the SU{2) x f/(l)-symmetry by the higgs vacuum condensate (vacuum expectation 
value - vev) t] and the W and Z bosons became massive. Let us consider the bare coupling 
constants (the bare charges cq of the photon, go of the PF-boson and /o of the Z-boson) and 
the bare masses of the vector bosons: 

mzo = ^foV , (75) 

mwo = ^9oV , (76) 

and also bare masses: rrito of the t-quark and mno of the higgs. 

Step II. We express a, G^, mz in terms of /o, 9o, Cq, tj, nito, niHo and 1/e (see Appendix D). 
Here 1/e appears because we use the dimensional regularization, calculating Feynman integrals 
in the space of D dimensions (see Appendix A). These integrals diverge ai D = 4 and are finite 
in the vicinity of D = 4. By definition, 

2e = 4:-D ^0 . (77) 

Note that in the one- loop approximation rritQ = rrit, mno = fnfj, since we neglect the electroweak 
corrections to the masses of the t-quark and the higgs. 

Step II is almost physics: we calculate Feynman diagrams (we say 'almost' to emphasize that 
observables are expressed in terms of nonobservable, 'bare', and generally infinite quantities). 

Step III. Let us invert the expressions derived at step II and write /o, 9o, V in terms of a, 
G^, TJiz, rrit, rriH and 1/e. This step is a pure algebra. 

Step IV. Let us express Vm, Va, Vr (or the electroweak one-loop correction to any other 
electroweak observable, all of them being treated on an equal basis) in terms of /o, 9o, r], rrit, 
rriH and 1/e. (Like step II, this step is again almost physics.) 

Step V. Let us express Vm, Va, Vr (or any other electroweak correction) in terms of a, 
G^, mz, m,t, m,H using the results of steps HI and IV. Formally this is pure algebra, but in 
fact pure physics, since now we have expressed certain physical observables in terms of other 
observables. If no errors were made on the way, the terms 1/e cancel out. As a result, we arrive 
at formula ( |52D which gives Vi as elementary functions of t, h and s. 

The five steps outlined above are very simple and visually clear. We obtain the main re- 
lations without using the 'heavy artillery' of quantum field theory with its counterterms in 
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the Lagrangian and the renormahzation procedure. This siniphcity and visual clarity became 
possible owing to the one-loop electroweak approximation (even though this approach to renor- 
mahzation is possible in multiloop calculations, it becomes more cumbersome than the standard 
procedures). As for the QCD-corrections to quark electroweak loops hidden in the terms (5VJ 
in equation (|52D, we take the relevant formulas from the calculations of other authors. 



5 One-loop corrections to hadronic decays of the Z bo- 
son. 

5.1 The leading quarks and hadrons. 



As discussed above (see formulas (pl|)- (|37|) ), an analysis of hardronic decays reduces to the 
calculation of decays to pairs of quarks: Z — > qq. The key role is played by the concept of 
leading hadrons that carry away the predominant part of the energy. For example, the Z ^r cc 
decay mostly produces two hadron jets flying in opposite directions, in one of which the leading 
hadron is the one containing the c-quark, for example, D~ = cd, and in the other the hadron 
with the c-quark, for example, Z)° = cu or A^ = udc. Likewise, Z ^ bb decays are identified 
by the presence of high-energy B ot B mesons. If one selects only particles with energy close 
to mz/2, the identification of the initial quark channels is unambiguous. The total number of 
such cases will, however, be small. If one takes into account as a signal less energetic 5-mesons, 
one faces the problem of their origin. Indeed, a pair bb can be created not only directly by a 
Z boson but also by a virtual gluon in, say, a Z —>■ cc decay ot Z —>■ uu, or ss. This example 
shows the sort of difficulty encountered by experimentalists trying to identify a specific quark- 
antiquark channel. Furthermore, owing to such secondary pairs, the total hadron width is not 
strictly equal to the sum of partial quark widths. 

We remind the reader that for the partial width Fg of the Z ^ qq decay we had eq. (|3l|) , 
where the standard width Fq was given by eq.(^) and the radiators RAg and Ryq are given in 
Appendix E. As for the electroweak corrections, they are included in the coefficients g^q and 
gvq- The sum of the Born and one- loop terms has the form 

9Aq = T,q[l + ^^^VAq{t,h)] , (78) 

310 I 

Rq = 9Vq/9Aq = 1 - MQ,\s^ + 4^ / J2 _ ^2) ^^^^ (^^ ^) " (^9) 

5.2 Decays to pairs of light quarks. 

Here, as in the case of hadronless observables, the quantities V that characterize corrections 
are normalized in the standard way: V —>■ t as t ^ 1. Naturally, those terms in V that are due 
to the self-energies of vector bosons are identical for both leptons and quarks. The deviation 
of the differences VAq — Vai and V^q — Vm from zero are caused by the differences in radiative 
corrections to vertices Z ^ qq and Z —^ II. For four light quarks we have 

1287rs^c^ 
VAu{t, h) = VAc{t, h) = VAi{t, h) + [—— {Fm + Fau) = 0.2634] , (80) 
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VAdit, h) = VAsit, h) = VAiit, h) + [—— {Fm - FAd) = 0.6295] 

6a 



VRu{t, h) = VUt, h) = Vm{t. h) + [ )- 1 X 

6a 

x[Fyi - (1 - As^)Fai + ^(-(1 - ^s^)Fau + Fvu)] = 0.1220] , 



^2) 



1 fiTTsrfr — s ") 

VRd{t, h) = Vnsit, h) = Vniit, h) + [ )- '- X 

6a 

x[Fvi - (1 - As^)Fai + 3((1 - P^)FAd - Fvd)] = 0.2679] , 



^3) 



where (see 41 



Fai = — (3.0099 + 16.45s2) , 

An 

Fvi = -J^(3.1878 + 14.9(5s^) , 

a 
Fau = - — (2.6802 + U.76s^) 
An 



a 



Vu 



-^(2.7329 + 14.2(5^2) 

An 



a 



FAd = ^{2.2221 + I6.^5s^) , 
An 



a 



Fvd = — (2.2287 +13.5(5s^ 
An 



U) 
^5) 
^6) 
^7) 
^8) 
^9) 



The values of F are given here for s^ = 0.23110 — Ss"^. The accuracy to five decimal places 
is purely arithmetic. The physical uncertainties introduced by neglecting higher-order loops 
manifest themselves already in the third decimal place. 

In addition to the changes given by eqs. (pOD - (p3D, one has to take into account also emission 
of virtual or "free" gluon from a vertex quark triangle. 

The corresponding effect cannot be parametrized in terms Vaq and V^g, because it con- 
tributes also to the radiators Raq and Rvq- The change of F/j caused by it has been calculated 
only recently |^6| and turned out to be rather small: 



STh{Z -^ u, d, s, c) = -0.59(3)MeV 



(90) 



5.3 Decays to 66 pair. 

In the Z -^ bb decay it is necessary to take into account additional t-dependent vertex correc- 
tions: 



VAb{t, h) = VAd{t, h) - 

VRb{t,h) = Vm{t,h)- 



is^c^ 



3^3 _ 2^2) (<^W + ^">W) = 5-03 



As\c^ - s') 
3(3 - 2s2) 



{4>{t) + 6^Mt)) = 1.76 



(91) 
(92) 
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Here the term 0(t) calculated in [^ corresponds to a ttW vertex triangle (see Figure 11 (a)), 



while the term Sag4>{t) calculated in |^, corresponds to the leading gluon corrections to the 



term 0(t) (see Figure 11 (b)): 6a^(j){t) ~ agt. Expressions for 0(t) and 6a^4'{t) are given in refs. 
511, m. For mt = 175 GeV, d,(m^) = 0.120 



0(t) = 29.96 , (93) 

5,Xt) = -3.02 , (94) 

and correction terms in equations (PTI) and (P2I) are very large. The subleading gluon corrections 



to (j){t) calculated recently |^ are very small: SThlZ ^ b) = —0.04 MeV. 

6 Comparison of one-electroweak-loop results and ex- 
perimental LEP-I and SLC data. 

6.1 LEPTOP code. 

A number of computer programs (codes) were written for comparing high-precision data of 
LEP-I and SLC. The best known of these programs in Europe is ZFITTER [Q, which takes 



into account not only electroweak radiative corrections but also all purely electromagnetic ones, 
including, among others, the emission of photons by colliding electrons and positrons. Some 
of the first publications in which the t quark mass was predicted on the basis of precision 
measurements |51]], were based on the code ZFITTER. Other European codes, BHM, WOH 



52 , TOPAZO [p3l, somewhat differ from ZFITTER. The best known in the USA are the 



results generated by the code used by Erler and Langacker IQ, [p^ . 

The original idea of the authors of this review in 1991-1993 was to derive simple analytical 
formulas for electroweak radiative corrections, which would make it possible to predict the t- 
quark mass using no computer codes, just by analyzing experimental data on a sheet of paper. 
Alas, the diversity of hadron decays of Z bosons, depending on the constants of strong gluon 
interaction a^, was such that it was necessary to convert analytical formulas into a computer 



program which we jokingly dubbed LEPTOP |5^. The LEPTOP calculates the electroweak 



observables in the framework of the Minimal Standard Model and fits experimental data so as 
to determine the quantities mt, mu and ds(wi^). The logical structure of LEPTOP is clear from 
the preceding sections of this review and is shown in the Flowchart. The code of LEPTOP can 
be downloaded from the Internet home page: [http://cppm.in2p3.fr./leptop/intro_leptop.htm] 



A comparison of the codes ZFITTER, BHM, WOH, TOPAZO and LEPTOP carried out in 



1994-95 [20| has demonstrated that their predictions for all electroweak observables coincide 



with accuracy that is much better than the accuracy of the experiment. The Flowcharts of 



LEPTOP and ZFITTER are compared on pages 25 and 27 of |20[; numerical comparison of 
five codes (their 1995 versions) for twelve observables is presented in figures 11-23 of the same 
reference. The results of processing the experimental data using LEPTOP are shown below. 

6.2 One-loop general fit. 

Second column of Table 4 shows experimental values of the electroweak observables, obtained 
by averaging the results of four LEP detectors (part a), and also SLC data (part b) and the 
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data on W boson mass (part c). (The data on the W boson mass from the pp-colhders and 
LEP-II are also shown, for the reader's convenience, in the form of s^r, while the data on s^ 
from i/A^-experiments are also shown in the form of mw- These two numbers are given in 
italics, emphasizing that they are not independent experimental data. The same refers to sf 
{^lr}-) We take experimental data from the paper |2^. The experimental data of Table 4 are 



used for determining (fitting) the parameters of the standard model in one-electroweak-loop 
approximation: rrit, rnn, o^sii^z) and a. (In fitting rrit the direct measurements of rrit by CDF 
and DO [collaborations] [^ are also used. In fitting a its value from eq.(|T3p was used.) The 
third column shows the results of the fit of electroweak observables with one loop electroweak 
formulas. The last column shows the value of the 'pull'. By definition, the pull is the difference 
between the experimental and the theoretical values divided by experimental uncertainty. The 
pull values show that for most observables the discrepancy is less than la. The number of 
degrees of freedom is 18 — 4 = 14. 
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Table 4. Fit of the experimental data [ J22[ | with one-electroweak-loop formulas 



mz = 91.1867(21) GeV is used as an input. 
Output of the fit: nin = ISQ.ltrel^^ GeV, d, = 0.1195 ± 0.0030, xV^rf.o./. = 15.1/14 





Experimental 


Fit 




Observable 


data 


standard 
model 


Pull 


a) LEP 








shape of Z-peak 








and lepton asymmetries: 








Tz [GeV] 


2.4939(24) 


2.4959(18) 


-0.8 


Ch [nb] 


41.491(58) 


41.472(16) 


0.3 


Ri 


20.765(26) 


20.747(20) 


0.7 


A'- 


0.0168(10) 


0.0161(3) 


0.8 


r-polarization: 








Ar 


0.1431(45) 


0.1465(14) 


-0.8 


Ae 


0.1479(51) 


0.1465(14) 


0.3 


Results for b and c 








quarks: 








Rfj* 


0.2166(7) 


0.2158(2) 


1 


Re* 


0.1735(44) 


0.1723(1) 


0.3 




0.0990(21) 


0.1027(10) 


-1.8 


^FB 


0.0709(44) 


0.0734(8) 


-0.6 


Charge asymmetry for pairs 








of light quarks qq: 








sKQfb) 


0.2321(10) 


0.2316(2) 


0.5 


b) SLC 








Alr 


0.1504(23) 


0.1465(14) 


1.7 


sKAlr) 


0.2311(3) 


0.2316(2) 


-1.7 


Rb* 


0.2166(7) 


0.2158(2) 


0.9 


Re 


0.1735(44) 


0.1723(1) 


0.3 


A, 


0.8670(350) 


0.9348(1) 


-1.9 


Ae 


0.6470(400) 


0.6676(6) 


-0.5 


c) pp+ LEP-11 +pN 








mw [GeV] {pp) + LEP-11 


80.39(6) 
0.2228(13) 


80.36(3) 


0.5 


s'w i^N) 


0.2254(21) 
80.2547(1089) 


0.2234(6) 


0.9 


mt [GeV] 


173.8(5.0) 


171.6(4.9) 


0.4 



Experimental values of Rb and Re correspond to the average of LEP-I and SLC results. 
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Table 5 



Observable 


sf 


Average over 
groups of observations 


Cumulative X^l^d.o.f. 
average 


ApB 

Ar 

Ae 


0.23117(55) 
0.23202(57) 
0.23141(65) 


0.23153(34) 


0.23153(34) 1.2/2 


^FB 

A^ 
^FB 


0.23226(38) 
0.23223(112) 


0.23226(36) 


0.23187(25) 3.4/4 


<Qfb > 


0.23210(100) 


0.23210(100) 


0.23189(24) 3.4/5 


Alr (SLD) 


0.23109(30) 


0.23109(30) 


0.23157(19) 7.8/6 



Table 5 [] gives experimental values of sf. The third column was obtained by averaging of 
the second column, and the fourth by cumulative averaging of the third; it also lists the values 
of x^ over the number of degrees of freedom. 



7 Two-loop electroweak corrections and theoretical un- 
certainties. 



In this Section we will discuss heavy top corrections of the second order in aw to mw and to 
coupling constants of Z-boson with fermions. Full calculation of a^r corrections is still absent. 

mt/mzY 

m. 



What have been calculated are corrections of the order ct^ t^ = a 
corrections ~ aLt 



56, 57 and 



vy- 



There are two sources of a^t^ corrections in our approach. The first source are reducible 
diagrams with top quark in each loop. The second source are irreducible two-loop Feynman 
diagrams which contain top quark |^, ^. We start our consideration with the first source the 
contribution of which is proportional to (11^(0) — Tlw{'^)Y ■ Detailed calculations are presented 
in Appendix F. 



7.1 OL^t^ corrections to mw /f^z , Qa and gv /qa from reducible 
diagrams. 



We start our consideration from the ratio of vector boson masses. From eq.( |l''.12 ) and ( |1''.13 ) 
we obtain: 



mw 
mz 



c[l + 



2(c2 - s^ 



, 3c^ - 10cS2 .^ 



Kc2 - S 



2^3 



Substituting the expression for 5 from ([F.10|) and using definition of Vm eq 



3(3-10s2)at2 



obtain the following correction to the function V„i'- 

47rsV(3-10s2) 

Oa V^, = --: r--^ -- 



4 "^ o./.2_ ^2p - 547^^2^^2-82)2 ' 



3a(c2 



(95) 
2D we 

(96) 



^Table 5 is our recalculation with LEPTOP program of the Table 30 of EWWG report [g2|. The numbers 
for Ae and Ar in tables 4 and 5 agree with each other, while they disagree in EWWG report in tables 30 and 
31. In order to restore the agreement one has to interchange Af. and Ar in table 30 in EWWG report. 
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The correction to axial coupling constant qai is easily derived from equations ( [b'.14| ) (since 
dAi ~ /o), ([F.10|) and definition of Vai, equations (|4^), (|5^): 



9Ai = -- - -S - —6^ , (97) 

Finally, taking into account definition of Vr, equations (|^), (|52D, and equations ( [F.15| ), 
we get: 

^y^Mz = 1 - 4[1 - c - -^^^s + j^zr^i^ ] = 

= 1 - 4^' + 4^^ - 7T^^' ' (99) 

c — s (c — S )'^ 

Formulas (p6D, (p8D and ( [L00| ) contain corrections to the functions \^ which come from the 
squares of polarization operators and are proportional to dt"^ - so, it is leading (~ t^) parts of 
{Hz — Hvi/)^ corrections. Numerically they are several times smaller than at'^ corrections which 
originate from irreducible diagrams. 

7.2 cx^t^ corrections from irreducible diagrams. 

The major part of a^t^ corrections comes from the irreducible two-loop Feynman diagrams 



56| , |57| . The key observation in performing their calculation is that these corrections are of 
the order of [c^wi^)^]'^ ~ ^t^ where Xt is the coupling constant of higgs doublet with the 
top quark. That is why they can be calculated in a theory without vector bosons, taking into 
account only top-higgs interactions [^ . Corresponding pieces of vector boson self-energies can 
be extracted from the self-energies of would-be-goldstone bosons which enter Higgs doublet 
(those components which after mixing with massless vector bosons form massive W- and Z- 
bosons). Correction of the order of a^t^ is contained in the difference nz(0) — n^(0) (see 
Figure 12), so it is universal, i.e. one and the same for Vm, Va and Vr. In [^ we call these 
corrections S^Vi: 

6Mt, h) = -T7-^^A{h/t) ■ f , (101) 

where function A{h/t) is given in the Table 6. To obtain this Table for rriH/fnt < 4 we use a 
Table from the paper |^, and for rriH/rnt > 4 we use expansion over rrit/mH from the paper 
m. For mt = 175 GeV and m^ = 150 GeV we get A = 6.4 and 64Vi{t,h) = -0.08. This 
corresponds to the shifts: —12 MeV for rriw, 7 ■ 10~^ for sf and 5 ■ 10~^ for qai- One should 
compare these shifts with one-loop results: ^iioop^^vK = 400 MeV, Snoopsf = 50 ■ 10~^ and 
^iioopfl'yi = 100 ■ 10^^. Let us remind that present experimental accuracy in mw is 64 MeV, in 
sf is 20 ■ 10^5 and in qai is 30 ■ 10~^ 

There is one more place from which corrections ~ a^t^ appear: this is the Z ^ bb decay. 
At one electroweak loop t-quark can propagate in the vertex triangle (tiW) (see Section 5). 
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That is why at two loops correction of the order a^t^ emerges. Due to this correction functions 
VAbit, h) and VRfe(t, h) differ from the corresponding functions describing Z ^ dd decay: 



^s^c^ 



VAbit, h) = VAd{t, h) - 3(3 _ 2,2) ('^(^) + ^'^(^' ^)) ' (102) 

VRbit, h) = V^t, h) - ^i^-Z_LZ(0(t) + 50(t, H)) , (103) 

where function 0(t) was discussed in Section 5 and 

50(t, h) = 5^M<P + SH<P{t, h) = ^ {-y (^)^ + ^^O'^'i^'^ ' ^^°^) 

First term in curly braces, 5o^0, was taken into account earlier, see Section 5, and the new 
correction 5H(p{t, h) is proportional to function r^ (h/t). Function r^ is given in Table 6. To 



obtain this Table for rriH/Tnt < 4 we use a Table from the paper |5^, and for mH/mt > 4 we 



use expansion over mt/niH from the paper [TO in full analogy with function A{h/t). 
For mt = 175 GeV, rriH = 150 GeV we have r^ = 1.6. 



The change of Tj, due to r^ = 1.6 equals 0.03 MeV, which corresponds to 2- 10 ^ shift in i?6. 



while experimental accuracy in i?^ is 7 ■ 10 (the one e-w loop correction in Ri, is —3.9 ■ 10 
The influence of r^ on A^p^ and Ai, is even smaller (by a few orders of magnitude) 

Table 6: Functions A{mH/'mt) and r^ {mn/'mt)- 
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rriH/mt 


A{mH/mt) 


(2) 

n '(mH/mt) 


rriH/mt 


A{mH/mt) 


(2) 

n '{mH/mt) 


.00 


.739 


5.710 


2.60 


10.358 


1.661 


.10 


1.821 


4.671 


2.70 


10.473 


1.730 


.20 


2.704 


3.901 


2.80 


10.581 


1.801 


.30 


3.462 


3.304 


2.90 


10.683 


1.875 


.40 


4.127 


2.834 


3.00 


10.777 


1.951 


.50 


4.720 


2.461 


3.10 


10.866 


2.029 


.60 


5.254 


2.163 


3.20 


10.949 


2.109 


.70 


5.737 


1.924 


3.30 


11.026 


2.190 


.80 


6.179 


1.735 


3.40 


11.098 


2.272 


.90 


6.583 


1.586 


3.50 


11.165 


2.356 


1.00 


6.956 


1.470 


3.60 


11.228 


2.441 


1.10 


7.299 


1.382 


3.70 


11.286 


2.526 


1.20 


7.617 


1.317 


3.80 


11.340 


2.613 


1.30 


7.912 


1.272 


3.90 


11.390 


2.700 


1.40 


8.186 


1.245 


4.00 


11.396 


2.788 


1.50 


8.441 


1.232 


4.10 


11.442 


2.921 


1.60 


8.679 


1.232 


4.20 


11.484 


3.007 


1.70 


8.902 


1.243 


4.30 


11.523 


3.094 


1.80 


9.109 


1.264 


4.40 


11.558 


3.181 


1.90 


9.303 


1.293 


4.50 


11.590 


3.268 


2.00 


9.485 


1.330 


4.60 


11.618 


3.356 


2.10 


9.655 


1.373 


4.70 


11.644 


3.445 


2.20 


9.815 


1.421 


4.80 


11.667 


3.533 


2.30 


9.964 


1.475 


4.90 


11.687 


3.622 


2.40 


10.104 


1.533 


5.00 


11.704 


3.710 


2.50 


10.235 


1.595 









7.3 CKwt corrections and the two-loop fit of experimental data. 



Corrections of the order aLt originate from the top loop contribution to W- and Z-boson self- 



2\2 



energies with higgs or vector boson propagating inside the loop and are of the order of g X 
We take into account these corrections in our code LEPTOP using results of the papers 

m. 
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Before we will present results of electroweak precision data fit which take into account a^ 
corrections, described in this Section, we must discuss how good the approximation which takes 
into account a^t^ and a^t terms but neglects (still not calculated) a^ terms should be. For 
rrit = 175 GeV we obtain t ^ 3.7, thus at first glance we have good expansion parameter so 
that a'yy terms could be safely neglected. To check this let us consider first the one electroweak 
loop, where the enhanced awt terms can be compared with non-enhanced aw terms. 

By using eqs.(^) and ( P9| ) and by comparing them with experimental data one sees that 
for mw /tTiz the awi term is equal to 0.0057, while the aw term is —0.0014. As for gvi/dAi, 
the two terms are 0.0122 and —0.0142. Thus for mw/^^z the awt term dominates, while for 
gvi/dAi it is practically cancelled by the aw term. 

Coming back to two-loop corrections we observe, that a^t^ correction to mw is not larger 
than a^^t correction; for rrit = 175 GeV and mn = 150 GeV it diminishes mw by 23 MeV 
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(compare with Section 7.2). 

In the Table 7 we present results of the fit of the data where we use theoretical formulas 
which include two-loop electroweak corrections described in this Section. Comparing Table 7 
with Table 4 where the fit of the one-loop electroweak corrected formulas was presented, we 
see that the fitted values of all physical observables are practically the same with one (very 
important) exception: the central value of the higgs mass becomes ~ 70 GeV lower. In view of 
the previous discussion it seems reasonable to consider this shift as a cautious estimate of the 
theoretical uncertainty in mn- 

We have a simple qualitative explanation why a'^t corrections reduce the higgs mass by 
~ 70 GeV. The point is that these corrections shift the theoretical value of sf by +0.0002, 
which is close to experimental error in sf. In order to compensate the shift the fitted mass of 
the higgs changes. This change can be easily derived. Indeed, from eqs. (|49|) , (|52|) , (^Sj) we get: 

^^' = -TTT^I ^c^SHr = -0.00086(5/7^ , (105) 

lD7r(c^ — s^) 

while from Table 3 we see that changing mn from 150 GeV to 100 GeV gives SHji = +0.236 
and 6s^ = -0.0002. 

In Figure 13 the dependence of x^ oii the value of higgs mass is shown separately with 
and without inclusion of SLD data (Z-decays into heavy quark pairs are taken into account on 
both plots). When all existing data are taken into account we get central value of higgs mass 



iriH = 71 GeV which is twenty GeV below direct bound [^ from LEP-II search: mn > 95 
GeV. However, uncertainty in the value of tuh extracted from radiative corrections is quite 
large, thus there is no contradiction between these two numbers. 

At the end of this Section we would like to make two remarks demonstrating that one should 
not take too seriously the central values of mn extracted from the global fits. 

First, if one disregards the FNAL measurements of rrit, then one obtains from the fit: 

mt = 160.7+^:^ GeV , 

mn = 30.3^?^:^ GeV . 

Such value oirriH is by 1.5 standard deviations below the lower bound from direct searches of 
LEP-II. ( Note also that the fitted value of the top mass nit is substantially lower than measured 
at FNAL). 

Second, as it was stressed in ref. [^, the values of sf extracted from different observables 
lead to very different central values of mn- For example, from SLAG data on Alr it follows 
that rriH = 25 GeV with 90% confidence interval from 6 GeV to 100 GeV. Even smaller values 
of niH follow from LEP measurement of Ap^: mu = 4 GeV (0.2 GeV < niu < 95 GeV at 
90% C.L.). As for other asymmetries measured at LEP, they lead to much heavier higgs: from 
A^PB, for example, tuh = 370 GeV (100 GeV < run < 1400 GeV at 90% G.L.). That is why 
the average of all these values of tjih seems to be not very reliable. 
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Table 7. Fit of experimental data 



niz = 91.1867(21) GeV is used as an input. 
Output of tlie fit: tuh = 70.8l^^ GeV*, d, = 0.1194 ± 0.0029, xV^d, 



with two-electroweak-loop formulas. 

15.0/14 



o.f. 



Observable 


Experimental 


Fit Standard 


Pull 




data 


Model 




a) LEP-I 








shape of Z-peak and 








lepton asymmetries: 








Tz [GeV] 


2.4939(24) 


2.4960(18) 


-0.9 


ah [nb] 


41.491(58) 


41.472(16) 


0.3 


Ri 


20.765(26) 


20.746(20) 


0.7 


4' 


0.0168(10) 


0.0161(4) 


0.7 


r-polarization: 








Ar 


0.1431(45) 


0.1467(16) 


-0.8 


Ae 


0.1479(51) 


0.1467(16) 


0.2 


results for heavy quarks: 








Rb** 


0.2166(7) 


0.2158(2) 


1.0 




0.1735(44) 


0.1723(1) 


0.3 


4^ 

^FB 


0.0990(21) 


0.1028(12) 


-1.8 


4<^ 
^FB 


0.0709(44) 


0.0734(9) 


-0.6 


charge asymmetry for pairs of 








light quarks qq: 








sKQfb) 


0.2321(10) 


0.2316(2) 


0.5 


b) SLC 








si {A,r) 


0.2311(3) 


0.2316(2) 


-1.6 


Alr 


0.1504(23) 


0.1467(16) 


1.6 


Rb** 


0.2166(7) 


0.2158(2) 


0.9 


Rr 


0.1735(44) 


0.1723(1) 


0.3 


Ab 


0.8670(350) 


0.9348(2) 


-1.9 


Ac 


0.6470(400) 


0.6677(7) 


-0.5 


c) pp+ LEP-II +uN 








mw [GeV] (pp+ LEP-II) 


80.3902(64) 
0.2228(13) 


80.3659(34) 


0.4 


s'w {^N) 


0.2254(21) 
80.255(109) 


0.2233(7) 


1.0 


mt [GeV] 


173.8(5.0) 


170.8(4.9) 


0.6 



* The most optimistic errors on Mh are obtained in the fit including a{DH) ^ ~ 128.923(36) |3Q| and 
as{PDG) = 0.1178(23) from low energy data ||. Such a fit gives tuh = 93^^? GeV, m* = 171.3 ± 4.8 
GeV, as — 0.1184 ± 0.0018, x^ /rid.o.f. = 15.2/14. However the systematic errors due to the model assumptions 
used in the calculations of as{PDG) and a{DH) are not easy to estimate. That is why we prefer to use the 
result with less optimistic assumptions leading to bigger error in mu- 
** Experimental values of Rb and Re correspond to the average of LEP-I and SLC results. 
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As can be seen from the Table 8, the LEPTOP fit is very close to the ZFITTER fit p2 
and to the fit by Erler and Langacker ||2^. This indicates that theoretical uncertainties are 
very small, except for the non-calculated part of the corrections, that is common to all three 
programs. 



Table 8. Comparison of the LEPTOP fit with the ZFITTER fit |22[ and with the 

fit by Erler and Langacker ||23|| . 



Observable 


Experimental 


Fit LEPTOP 


Fit EWWG 


Fit Erler- 




data 




ZFITTER 


Langacker* 


a) LEP-I 










Mz [GeV] 


91.1867(21) 


91.1867 fix. 


91.1865 


91.1865(21) 


Tz [GeV] 


2.4939(24) 


2.4960(18) 


2.4958 


2.4957(17) 


ah [nb] 


41.491(58) 


41.472(16) 


41.473 


41.473(15) 


Ri 


20.765(26) 


20.746(20) 


20.748 


20.748(19) 


^FB 


0.0168(10) 


0.0161(4) 


0.01613 


0.0161(3) 


Ar 


0.1431(45) 


0.1467(16) 


0.1467 


0.1466(15) 


A, 


0.1479(51) 


0.1467(16) 


0.1467 


0.1466(13) 


Rb 


0.2166(7) 


0.2158(2) 


0.2159 


0.2158(2) 


Re 


0.1735(44) 


0.1723(1) 


0.1722 


0.1723(1) 


ApB 


0.0990(21) 


0.1028(12) 


0.1028 


0.1028(10) 


4c 


0.0709(44) 


0.0734(9) 


0.0734 


0.0734(8) 


sKQfb) 


0.2321(10) 


0.2316(2) 


0.23157 


0.2316(2) 


b) SLC 










sKAlr) 


0.2311(3) 


0.2316(2) 


0.23157 


- 


Alr 


0.1504(23) 


0.1467(16) 


— 


0.1466(15) 


Ak 


0.8670(350) 


0.9348(2) 


0.935 


0.9347(1) 


A, 


0.6470(400) 


0.6677(7) 


0.668 


0.6676(6) 


c) pp+ LEP-II +vN 










mw [GeV] {pp+ LEP-II) 


80.3902(64) 

0.2228(13) 


80.3659(34) 


80.37 


80.362(23) 


s'w i^N) 


0.2254(21) 
80.255(109) 


0.2233(7) 


0.2232 




mt [GeV] 


173.8(5.0) 


170.8(4.9) 


171.1(4.9) 


171.4(4.8) 


mu [GeV] 




71 +82 
'-■-•-43 


76.+«^ 


107.1^^ 


tts 




0.1194(29) 


0.119(3) 


0.1206(30) 


a-i 


128.878(90) 


128.875 


128.878 





* Erler-Langacker use slightly different experimental datasct for their fit. This may cause some of the 
discrepancies with LEPTOP and ZFITTER. 
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8 Extensions of the Standard Model. 

The Standard Model works well at the energy scale of the order of the vector bosons masses. 
We see that the SM description of the electroweak observables in this energy region is in perfect 
agreement with the precision measurements. 

However there are many natural physical questions that have no satisfactory answers within 
the framework of the SM. So it is hard to believe that the Standard Model is the Final Theory. 
The common expectation is that there should be New Physics beyond the Standard Model. 

Direct accelerator searches did not find yet any trace of New Physics. Their negative results 
gave lower bounds on the masses and upper bounds on the production cross sections for the 
new particles. In this section we are going to study the indirect bounds on New Physics that 
can be theoretically derived from the precision measurements at low energy of the order of Z 
and W boson masses. Loops with hj^othetical new particles change the predictions of the SM 
for electroweak observables. Since the SM gives a very good description of the data there is 
little room for such new contributions. In this way one can get some kind of constraints on 
new theory. 

Any possible generalizations of the SM are naturally divided into two classes: theories with 
and without decoupling. In the first class the contribution of new particles into W and Z boson 
parameters are suppressed as positive powers of {rri^/m'^)'^ when the masses of new particles 
m become larger than electroweak scale. One cannot exclude such theory by studying loop 
corrections to low-energy observables. In this way one may have hopes to bound the masses 
of new particles from below. The most famous example of such theory are supersymmetric 
extensions of the SM. 

In the second class of theories the contribution of new particles into low-energy observables 
does not decouple even when their masses become very large. Such SM generalizations can be 
excluded if the additional nondecoupled contributions exceed the discrepancy between SM fit 
and experimental data. The example of such generalization is the SM with additional sequential 
generations of quarks and leptons. 

8.1 Sequential heavy generations in the Standard Model. 

We start the discussion of New Physics with the simplest extension of the SM, namely with the 
SM with additional sequential generations of leptons and quarks ([^- [H])- Nobody knows 
any deep reason for the number of generations to be equal to three. So it is interesting to 
study whether it is allowed to have four and more generations. Certainly these new generations 
should be heavy enough not to be produced in the Z decays and at LEP-II. 

We consider the case of no mixing between the known generations and the new ones. In 
this case the new fermion generations affect the ratio itiw/ttiz and the widths and the decay 
asymmetries of the Z boson only through the vector bosons self-energies. Such kind of correc- 
tions have been dubbed [|^ "obhque corrections". We start their study with the case of SU{2) 
degenerate fourth generation: 

mu = rriD = rriQ , tun = tue = tul (106) 

New terms in the self-energies modify the functions Kn, Va, Vr , i.e. the radiative corrections 
to mw /tuz, Qai and gvi/dAi- The contribution to Vi from the fourth generation can be written 



32 



in the form: 

Vm^Vra + 5^V^ , Va^Va + 5^Va , Vr^Vr + 5^Vr . (107) 

The analytical expressions for d'^Vi for quark or lepton doublets (neglecting gluonic corrections) 
can be found in Appendix G, eqs (G.1)-(G.3). 

In the limit of a very heavy fourth generation of leptons and quarks one has : 

^5'Vm^-^s\ S5V«^-^, S^Va^O , (108) 

where S denotes sum over leptons and quarks with mq = ttil = 1714^, s^ ~ 0.23. Equations 
( |108| ) reflect the non-decoupling of the heavy degrees of freedom in electroweak theory, caused 
by the axial current. It is interesting that the contribution of degenerate generation to Kn, Vr 
has negative sign. 

The fourth generation with strong violation of SU{2) symmetry (i.e. with very large mass 
difference in the doublet) gives universal contribution to functions 5^Vi (similar to the universal 
contribution of t- and 6-quarks from the third generation to Vj) 

S^Vi = \\m% - m||/3m|. (109) 

In the case of large mass splitting (5^V^ are positive. From eqs. ( |108| ) and ( |109| ) it is clear that 
somewhere in the intermediate region of mass splitting the functions (5^Kn and ^^Vr intersect 
zero. In the vicinities of these zeroes the contribution of new generation to these specific 
observables is negligible and one can not exclude these regions of masses studying only one of 
the observables. Fortunately for different observables these zeroes are located in the different 
places and the general fit overcomes such conspiracy of new physics. 

For different up- and down- quark (and lepton) masses analytical expressions for 8'^Vi are 
given in Appendix G, eqs. (|G.4| ) - (|G.6|) . 



Figure 14 demonstrates 2-dimensional exclusion plot for the case of n extra generations, 
where n is formally considered as a continuous parameter. We see from this plot that at 90% 
c.l. we have less than one extra generation and at 99% c.l.- less than two extra generations for 
any differences of up- and down- quarks masses. 

8.2 SUSY extensions of the Standard Model. 

In this section we consider another example of new physics: supersymmetric extensions of the 
SM. There are certain aesthetic and conceptual merits of such SUSY generalization of the SM. 
Here are some of them: 

1) Supersymmetry gives a solution for the problem of fine tuning, i.e. it prevents the 
electroweak scale of the SM from mixing with the Planck scale. 

2) The problem of unification of electroweak and strong coupling constants seems to have 
solution in the framework of SUSY extensions. 

3) Finally, any ambitious "Theory of Everything" inevitably includes SUSY as the basic 
element of the construction. 

To make systematic introduction into SUSY extensions of the SM one needs a separate 
review paper ( for the review papers see e.g. [^ ). Here we are going to make a short 
sketch of this well developed branch of physics in applications to the theory of Z boson. To 
construct SUSY extensions one has to introduce a lot of new particles. For example minimal 
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N = 1 supersymmetry automatically doubles the number of degrees of freedom of the SM: any 
fermionic degree of freedom has to be coupled with bosonic degree of freedom and vice versa. 
Thus the left ( right ) leptons have to be accompanied by scalars: "left" ("right" ) sleptons, 
quarks by squarks, gauge bosons by spinor particles - gauginos, etc. The Higgs mechanism 
of mass generation for up and down quarks requests for two Higgs boson doublets ( and two 
higgsino doublets respectively). 

Not a one of these numerous new particles has been observed yet. If they do exist they 
are too heavy to be produced at the working accelerators. On the other hand, these heavy 
supersymmetric particles (again if they do exist) are produced in the virtual states, i.e. in 
the loops. Loops with new particles change the predictions of the SM for the low-energy 
observables. ( Under "low-energy" we mean here E < mz )■ In this indirect way one can get 
some information about existence or nonexistence of SUSY. 

The SUSY extensions of the SM belong to the class of new physics that decouples from 
the low-energy observables when the mass scale of this new physics becomes very large. It 
means that the additional contribution into electroweak observables due to the supersymmetric 
particles are of the order oi awifnw/^T^susY)"^ or avK('^t/''^5C/s'y)^ ; where msusY characterizes 
the mass scale of superpartners. Since the fit of the precision data in the framework of the SM 
statistically is very good these new additional contributions have to be small. So in this way 
one expects to get strong restrictions on the value of msusY- 



Supersymmetric contributions into low-energy observables were studied in papers ||6^ - |^ 
The results depend on the model and on the pattern of SUSY violation. Within a given model 
the results for low-energy observables are formulated in terms of the functions that depend on 
the fundamental parameters of the SUSY Lagrangian that are fixed at the high energy scale 
of SUSY violation. The fit of experimental data in the framework of a given SUSY model 
imposes certain restrictions on the allowed region of these high-energy scale parameters of the 
model. As for the masses of sparticles their values are calculated by numerical solution of the 
renormalization group equations. They also depend on the fundamental SUSY parameters at 
high energy scale. In this rather indirect way one gets restrictions on the physical masses of 
sparticles in general case. 

To give the reader the taste of exploration of the new supersymmetric physics we consider 
in this section only that part of the multi-dimensional space of SUSY parameters for which all 
sparticles have more or less the same masses, i.e. when we have no light sparticles. (It seems 
reasonable to start the study of the unknown field with such kind of the simplest assumptions). 
In this case one can find the class of enhanced oblique corrections that are universal, i.e. that 
are the same for any model. Another merit of these corrections is that they directly depend on 
the masses of sparticles. 

As will be shown, the enhanced electroweak radiative SUSY corrections are induced by the 
large violation of SU{2)l symmetry in the third generation of squarks. Therefore we start the 
discussion of the SUSY corrections to the functions Vi with the brief description of the stop 
{tL,tR) and sbottom {pLihR) sector of the theory. The following relation between masses of 
quarks q and diagonal masses of left squarks qi takes place in a wide class of SUSY models: 

ml^ =ml + mlusY + "^1 cos(2/5)(r3 - s^Q,) , (110) 

where s^ ^ 0.23, Qq is the charge and T3 is the third projection of weak isospin of quark and tg/3 
is equal to the ratio of the vem of two Higgs fields, introduced in SUSY models. The second 
term in the r.h.s of the eq.( |110D violates supersymmetry. It is some universal SU(2)-blind 



34 



SUSY violating soft mass term. The third term in r.h.s of the eq.( |llU[ ) also violates SUSY. 
It originates from quartic D term in the effective potential and is different for up and down 
components of the doublets. The only hypothesis that is behind this relation is that the origin 
of the large breaking of this SU{2)l is in the quark-higgs interaction. 

Therefore from eq. (|110|) we get the following relation between masses of stop m? , of sbottom 



rrij and of top mf ( we neglect nib 



1 1 

nir — rrii 



ml + m| cos(2(3)(? 



'1111 



Relation (|111|) is central for this approach. It demonstrates the large violation of SU{2)l 
symmetry in the third generation of squarks. On the other hand it demonstrates that in the 
limit of very large mass the left stop and left sbottom become degenerate and the parameter 
(m? —^ )/'^ goes to zero when m^f/sy goes to infinity. That is why the physical observables 
can depend on this decoupling parameter. 

As for the right sparticles from the third generation, they are SU{2)l singlets. But they 
can mix with the left sparticles and in this way they contribute into enhanced corrections. The 
mixing between Bl, 6r has to be proportional to nib and can be neglected. The tLtR mass matrix 
in general has the following form: 



7n} 



tL 

mtA'. 



mtA. 



niz 



(112) 



where titR mixing is proportional to nit and therefore is not small. Coefficient A'f- depends on 



the model. Diagonalizing matrix ( p.l2| ) we get the following eigenstates: 



^1 — CutL + SutR^ 
^2 = ~SutL + CutR 



(113) 



where Cu = cos 9lr, Su = sin 9lr, Olr is the titR mixing angle, and 



tg^^ 



m^ 



m 



LR 



tL 



2 2 ' 

??lr — ni2 



2 \ 2 \ 2 

m^ > ni^ > ni2 ■ 



tL 



(114) 



Parameters mi and ni2 are the mass eigenvalues: 



m?2 



771? + m? 

tL tji 



± 



I 2 

\nii 

' tL 



TUz 



■\ 



Ani^A'i 



[m 



tL 



tR' 



(115) 



The enhanced electroweak radiative corrections are induced by the contribution of the third 
generation of squarks into self-energy operators of vector bosons. Nondiagonal vector currents 
of squarks are not conserved only because of violation of SU(2) by mass terms. Thus one should 
expect that the self-energy operators are proportional to the divergency of the currents. To 
calculate these enhanced terms it is sufficient to expand the operators of the vector bosons 
Sv'(fc^) at fc^ = 0. The terms enhanced as mf/Mg^gY come from Svi/(0), while those enhanced 
as mfM^^/Mg^gY come from E'^^(O) (see Figure 15). These simple self-energy corrections 
are obviously universal since stop and sbottom should exist in any SUSY model and the cou- 
pling constants are universal since they are fixed by gauge invariance only. The higher-order 



35 



derivatives of self energies are suppressed as {rnw,zl'fnsusYY- They are of the same order of 
magnitude as the numerous mo del- dependent terms coming from vertex and box diagrams. If 
there are no very light sparticles the first two universal terms have rather large enhancement 
factor of the order of t^ ~ 14 and t ^ 3.7 respectively. (The presence of terms ~ rnl in SUSY 
models was recognized long ago |7^ ). We neglect the non-enhanced terms. The accuracy of 
such approximation may be of the order of ten percent if we are lucky, but it may be as well of 
the order of unity (see discussion of Vr in Sections 5 and of the two-loop corrections in Section 
7). As for the stop contributions to the vertex corrections there is only one relevant case - the 
amplitude oi Z ^ bb decay. For vertex with stop exchange there are no terms enhanced as 
{mt/mwy [Ol- Thus we will neglect corresponding corrections as well. 

The calculation of the enhanced two terms is a rather trivial exercise. The only subtle 
point is the diagonalization of the stop propagators. The result of calculations depends on 3 
parameters: rrii, m2 and mi . The dependence on angle (3 is very moderate and in numerical 
fits we will use rather popular value tg/? = 2. In what follows instead of rrii we will write mi, 
bearing in mind that bibn mixing is proportional to m,h and can be neglected. The formulas 
that describe the enhanced SUSY corrections to the functions V^ can be found in the Appendix 
G, eqs.(^ - (IGTTl) . 



There is also another source of the potentially large SUSY corrections: vertices with gluino 
exchange of the order dsi^mz/m^susv)^- 



These corrections shift the radiators Rv„ and Ra„ in eq.(0) ||73| : 



SRv, = SRa, = 1 + ^^^^A,{x,y) , (116) 

TT 



1 1-21 



Ai(x,y) = -- / dzx \ rfzslogfl —. ■ TT r] , (117) 

V 1 Z/V3 / 

where x = (m^/mg)^, y = {mz/m.gY, and Ai(a;, x) ^ ■^x + .... We take these gluino corrections 
into account in our analysis. The electroweak SUSY corrections to qa and gy are generated 
by the corrections to the function Va eq.( p.7| ) and Vr eq.( |G.8| ). 



Having all the necessary formulas in hands we start the new fit of the data with the simplest 
case of the absence of tit^ mixing, syhOlr = 0. In this case we have only one additional mass 
parameter. Thus we expect that this mass should be heavy enough not to destroy the perfect 
SM fit of the experimental data. First let us take the value of the lightest neutral Higgs boson 
mass as a free parameter and take the masses of the other three Higgs bosons to be very heavy. 
The results of the fit are shown in Table 9. We see that to fit the data with light sbottom 
one has to take the mass of the Higgs boson much larger than its Standard Model fit value. 
Even in this case the quality of the fit is worse than the SM one. For very heavy sbottom one 
reproduces the SM fit. (To reduce the number of parameters we take m^ = m^ in this fit. Let 
us stress that light squarks with masses of the order of 100 — 200 GeV are usually allowed only if 
gluinos are heavy, m^ > 500 GeV [|7i|. In the case of heavy gluino the correction Ai (eq. ( |117| )) 
becomes power suppressed and we return to the Standard Model fit value of d^ = 0.119(3)). 
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Table 9 



m~, (GeV) 


m, (GeV) 


&s 


X^/nd.o.f. 


100 
150 
200 
300 
400 
1000 


850^ 
484+ii 

i52j:^f 

77+^^ 


0.113 ±0.003 
0.116 ±0.003 
0.117 ±0.003 
0.118 ±0.003 
0.119 ±0.003 
0.119 ±0.003 


20.3/14 
18.1/14 
17.3/14 
16.3/14 
15.8/14 
15.2/14 



Fit of the precision data with SUSY corrections taken into account in the case of the absence 
of tit n mixing, sin 6 lfi. = and nih taken as a free parameter. For m^ > 300 GeV SUSY 
corrections become negligible and SM fit of the data is reproduced. 

We see that to get reasonably good fit of the data in the framework of the SUSY extensions 
with light squarks one has to put the lightest Higgs mass in the TeV region. It is time to remind 
that in SUSY models the mass of Higgs boson is not an absolutely free parameter. In MSSM 
(Minimal Supersymmetric Standard Model) among three neutral Higgs bosons the lightest one 
should have mass less than approximately 120 — 135 GeV ||76|. If other higgses are considerably 
heavier the lightest scalar boson has the same couplings with gauge bosons as in the Standard 
Model. As a result the same SM formulas for radiative corrections can be used in the SUSY 
extensions of the SM. (Deviations from the SM formulas are suppressed as (m/i/m^)^, where 
niA - the mass of the heavier higgs. We will assume in our analysis that rriA is large). For 
maximal allowed value mh — 120 GeV the results of the fit are shown in Table 10. (In what 
follows we will always take uih — 120 GeV since for 90 GeV < rrih < 135 GeV results of the fit 
are practically the same.) This table demonstrates that superpartners should be heavy if we 
want to have good quality fit of the data. 



Table 10 



mi (GeV) 


&s 


xV^d.o./. 


100 


0.110 ±0.003 


30.2/15 


150 


0.115 ±0.003 


21.9/15 


200 


0.116 ±0.003 


18.6/15 


300 


0.118 ±0.003 


16.4/15 


400 


0.119 ±0.003 


15.8/15 


1000 


0.119 ±0.003 


15.5/15 



The same as Table 9 but with the value of the lightest Higgs boson mass mh 
which is about the maximum allowed value in the simplest SUSY models. 



120 GeV 



The next step is to take into account tit^ mixing. In the Figure 16 we show the dependence 
of SUSY corrections ^susvVi on mi and m2 for mi = 200 GeV. One clearly sees from this Figure 
that even for this small value of m^ there exist the domain of low m2 values where the enhanced 
radiative corrections are suppressed. In Figure 16 one sees the valley where SsusyVz reaches 
its minimum values which are considerably smaller than 1. The valley starts at m2 — m^, 
mi ~ 1000 GeV and goes to m2 — 100 GeV, mi ~ 400 GeV. The smallness of the radiative 



corrections at the point m2 



~ mi, nil 



1000 GeV can be easily understood: here 6lr ~ 7r/2, 



t2 — t^, ti c^ tji- Thus nondiagonal charged left current of squarks is conserved and the main 
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enhanced term vanishes. Indeed in SsusyVa only the term proportional to g{m2,mi) remains 



in eq.( |G.7|) , but for 1712 = m^ it is equal to zero. At this end point of the valley t2 —ti, ti ^ tR, 



, > rrij . Almost 
which means that m| > m~ . This possibility to 



so mj S> mf , which is opposite to the relation between mi^ and mi^ occurring in a wide class 

of models. In these models (e.g. in the MSSM) the left and the right squark masses are equal 

at high energy scale. When renormalizing them to low energies one gets m? ^ -^^ 

along the whole valley we have t^OiR > 1 

suppress radiative corrections was discussed in [[7S] . However, in the vicinity of the end point 

mi ~ 300 GeV, m2 — 70 GeV the value of tg^OiR becomes smaller than 1 and mj < mj . 

In Table 11 we show values of x^ along the valley of its minimum, which is formed for 
mi = 200 GeV. We observe that good quality of fit is possible for light superpartners if iitR 
mixing is taken into account. 



Table 11 



mi (GeV) 


ma (GeV) 


as 


xV^d.o./. 


1296 


193 


0.118 ±0.003 


15.6/15 


888 


167 


0.118 ±0.003 


15.8/15 


387 


131 


0.118 ±0.003 


16.1/15 


296 


72 


0.117 ±0.003 


16.7/15 



Results of fit along the valley of minimum of x^ for fixed value m^ ^ 200 GeV andrrih — 120 
GeV . 

The main lesson of this subsection is the following. The fit of the precision data on elec- 
troweak observables (i.e. of Z-boson decay parameters from LEP and SLC and the values of 
ruw and the rrit from Tevatron) in the framework of SUSY extension of the SM assuming small 
value of mi, the absence of tLtR mixing and m/i = 120 GeV leads to the growth of x^ value. For 
heavy squarks the SUSY sector of the theory decouples from low-energy observables and the 
results of Standard Model fit are reproduced. On the other hand even for light sbottom and for 
small mass of one of two stops, one can find the values of t^tR mixing where supersymmetric 
corrections appear to be small and not excluded by experimental data. In this case the quality 
of the fit (i.e. the value of x^) is almost the same as in the Standard Model. 

9 Conclusions. 

The comparison of LEP-I and SLC precision data on Z boson decays with calculations based 
on the Minimal Standard Model has confirmed the predictive power of the latter: 

1. It was proved that there exist only three generations of quarks and leptons with light 
neutrinos. 

2. The Z boson couplings of quarks, charged leptons and neutrinos are in accord with the 
theory. 

3. From the analysis of the radiative corrections the mass of the top quark had been correctly 
predicted before this particle was discovered at Tevatron. 
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4. All electroweak observables (except for the mass of the higgs) are perfectly fitted by one 
loop electroweak corrections (with virtual and "free" gluons being taken into account). 

5. The dependence of the radiative corrections on the mass of the higgs is feeble when higgs 
is heavy. Therefore the value of the higgs mass extracted from LEP-I and SLC data has 
rather large error bars. Within one standard deviation the central fitted value of the higgs 
mass becomes smaller when the leading two-electroweak-loop corrections are taken into 
account. In this case it is close to 90 GeV - its direct lower limit from the LEP-II search. 
However the non-leading two loop corrections may change this result. Calculation of these 
corrections is a challenge to theorists. Better understanding of systematic discrepancies 
between various asymmetries in Z-decays is a challenge to experimentalists. 

6. One of the main conclusions of the one-electroweak-loop case is that in this case the value 
of the leading and non-leading corrections are comparable and even may cancel each other 
(in the case of leptonic parity violating parameter gvil Qai)- 

7. The remarkable agreement between the Minimal Standard Model and experimental data 
on Z-decays puts strong limits on the hypothetical " new physics" , such as extra genera- 
tions of heavy quarks and leptons and/or properties of supersymmetric particles. 

The discovery of the higgs, more precise measurements of the mass of the W boson at LEP- 
II and Tevatron, and more accurate prediction of the value of electric charge at the scale of m^, 
may substantially improve the sensitivity of the Z-decay data to the possible manifestations of 
the new physics. 

We conclude this review with a flowchart summarizing the theoretical approach used by us. 
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Flowchart. 



Choose three observables measured with the highest accuracy: 

G^ ,mz , a{mz) = a 




Determine angle 9 ( s = sin^, c = cos 6*) 
in terms of G^ , mz , a: G^ = (7r/\/2)a/s^c^m| 




Introduce bare couphng constants in the framework of MSM 

{ao , azo , awo), bare masses (m^o, mwo, mno, rrito, m^o) 
and the vacuum expectation value (vev) of the higgs field, r/. 




Express ao , azo ? ^zo in terms of G^ , mz , a in one-loop 
approximation, using dimensional regularization (l/e ,ii). 




Express one-electroweak-loop corrections to all electroweak 

observables in terms of ao , azo , 'mzo , TUt , mn, and hence, in terms of 

of G^ , iJiz ^a^rrit, tuh- Check cancellation of the terms [l/e , fi). 




Introduce gluon corrections to quark loops and QED (and QCD) 

final state interactions, as well as the two-electroweak-loop 

corrections, calculated by other authors, in terms of 

a ,as{mz) ,mH ,mt , m^. 




Compare the predictions of the successive approximations 

(Born, one loop, two loops) with the experimental data 

on Z-decays. Perform the global fits for mn , rrit , as{mz) , a 

for one and two loops. Derive theoretical predictions 

of the central values for all electroweak observables 
and of the corresponding uncertainties ('errors'). 




Check the sensitivity of the fits to the possible 

existence of new heavy particles: extra generations, 

SUSY, technicolor, Z' etc. 
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Appendix A. 
Regular izat ion of Feynman integrals. 

Integrals corresponding to diagrams with loops formally diverge and thus need regulariza- 
tion. Note that there does not exist yet a consistent regularization of electroweak theory in all 
loops. A dimensional regularization can be used in the first several loops; this corresponds to 
a transition to a D-dimensional spacetime in which the following finite expression is assigned 
to the diverging integrals: 



vrT r(f + s)r(a- 


D 

2 


-)., 


r(f) r(a) 

(m2)f-"+« 






/i^-4 





where /i is a parameter with mass dimension, introduced to conserve the dimension of the 
original integral. 

This formula holds in the range of convergence of the integral. In the range of divergence, 
a formal expression ( |A.1| ) is interpreted as the analytical continuation. Obviously, the integral 



allows a shift in integration variable in the convergence range as well. Therefore, a shift p -^ p+q 
for arbitrary D can also be done in ( |A.1| ). This factor is decisive in proving the gauge invariance 



of dimensional regularization. 

At D = 4 the integrals in ( |A.1| ) contain a pole term 



2 m^ 

A = - - + log47r-7-log— , (A.2) 

where 7 = 0.577... is the Euler constant. Choice of constant terms in ( |A.2| ) is a matter of 
convention. 

The algebra of 7-matrices in the D-dimensional space is defined by the relations 

7/.7v + lulf, = '2^9 t,v X / , (A. 3) 

g,, = D , (A.4) 

7^7.7/. = (2 - ^)7. , (A.5) 

where / is the identity matrix. 

As for the dimensionality of spinors, different approaches can be chosen in the continuation 
to the /^-dimensional space. One possibility is to assume that the 7 matrices are 4x4 matrices, 
so that 

Spl = 4: . (A.6) 

The D-dimensional regularization creates difficulties when one has to define the absolutely 
antisymmetric tensor and (or) 75 matrix. For calculations in several first loops, a formal 
definition of 75, 

757/. + 7^75 = , (A. 7) 

ll = I (A.8) 
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does not lead to contradictions. 

Thus, the amphtudes of physical processes, once they are expressed in terms of bare charges 
and bare masses, contain pole terms ~ 1/(D — 4). 

If we eliminate bare quantities and express some physical observables in terms of other 
physical observables, then all pole terms cancel out. The general property of renormalizability 
guarantees this cancellation. (We have verified this cancellation directly in |^.) The "Five 
steps" described in Section 4.6 are based on this renormalization procedure. 

In order to avoid divergences in intermediate expressions, one can agree to subtract from 

each Feynman integral the pole terms ~ 1/(4 — D), since they will cancel out anyway in the final 

expressions. Depending on which constant terms (in addition to pole terms) are subtracted from 

the diagrams, different subtraction schemes arise: the MS scheme corresponds to subtracting 

the universal combination 

2 
J— ^-7 + log47r . 
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Appendix B. 
Relation between a and a(0). 

We begin with the following famous relation of quantum electrodynamics [^ : 

Here the fine structure constant a = q;(0) is a physical quantity. It can be measured as a 
residue of the Coulomb pole 1/g^ in the scattering amplitude of charged particles. As for the 
running coupling constant (y.{q^)-, it can be measured from the scattering of particles with large 
masses m at low momentum transfer: m ^ y\(f\- In the standard model we have the Z-boson, 
and the contribution of the photon cannot be identified unambiguously if g^ 7^ 0. Therefore, 
the definition of the running constant a{q^) becomes dependent on convention and on details 
of calculations. 

At q^ = m|, the contribution of VT-bosons to a = Oi{rn?z) i^ ^^ot large, so it is convenient to 
make use of the definition accepted in QED: 

where 

fe = -n^(m|) + e;(o) , 

n,(m|) = ^S,(m|) . (B.3) 

One-loop expression for the self-energy of the photon can be rewritten as |78|: 



S^(s) = (a/37r)^iV/Qj[sA/ + (s + 2m5)F(s,m/,m/)-s/3]- 
/ 
- (a/47r)[3sAvy + (3s + 4m^)F(s,mty,miy)] , (B.4) 

where s = q"^, the subscript / denotes fermions, the sum Sj runs through lepton and quark 
flavors, and N;[ is the number of colors. The contribution of fermions to S^(g^) is independent 
of gauge. The last term in ( p.4| ) refers to the gauge-dependent contribution of VF-bosons; the 
't Hooft-Feynman gauge was used in equation ( [B.4|) . 
The singular term Aj is: 

1 m- 

A, = --7 + log47r-log^ , (B.5) 

where 2e = 4 — D {D is the variable dimension of spacetime, e ^ 0), 7 = — r'(l) = 0.577... is 
the Euler constant and /i is an arbitrary parameter. Both 1/e and n vanish in relations between 
observables. 

The function F{s,mi,m2) is defined by the contribution to self-energy of a scalar particle 
at g^ = s, owing to a loop with two scalar particles (with masses mi and 1712) and with the 
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coupling constant equal to unity: 



772? + ml mi 

F{s,mi,m2) = -1 H ^ 2 ^^S — 

mf — m,2 m,2 



— / dxlog 




x'^s — x{s + m^ — m,l) + m^ — ie 
m,im,2 



:b.6) 



The function F is normalized in such a way that it vanishes at g^ = 0, which corresponds to 
subtracting the self-energy at g^ = 0: 



F(0,mi,m2) = , 
The following formula holds for m,i = m,2 = m,: 



F{s,m,m,) = F{t) 



1 - V4r - 1 arcsin -^ , 4r > 1 , 



l-v^r^47logi±^l , 4r<l , 



4r 



:b.7) 



(B.8) 



where r = m^ / s. 

Let us present the following useful equality which holds for F{t) derivative: 



(B.9) 



To calculate the contributions of light fermions, the t-quark and the H^-boson to 5a, we 
need the asymptotics F{t) for small and large r: 



F{t) ^logr + 2 
Fir) ' ' 



Qt 60r^ 



\r\ < 1 



+ -, kl>i , 



(B.IO) 
(B.ll) 



As a result we obtain 



n^(m^) 



S7(-l)_ "^Ar/g^(A^ + | 



772; 



a 



TT 



gj 



Sn 



At + (1 + 2t)F(t) - 



1 



where t = m^i/mr^, and 



- — [3Ah. + (3 + 4c^)F(c^)] 
47r 



S;(0) = ::^EA^cV/A/-^(3Ah^ + | 



37r 



A-n 



(B.12) 



(B.13) 



(5q; 



q; 

7r 



E 

. 8 



^^ log — 
3 V m 



5^ 
3 



- Ql 



[I + 2t)F{t) 



1 



3 
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+ 



+ c' Fie' 



(B.14) 



Therefore, 5a is found as a sum of four terms, 

6a = 6ai + 6ah + Sat + Saw 
In one-loop approximation: 



6ai = 
Higher loops |2^ give: 

Loops with top quarks give: 

5at ~ 



a 
3^Y 


ml 5 

log ^ - 7T 
mf S 


Sai = 0.031498 


a 4 
7r45 


V Mt ) 



0.03141 



0.00005(1) 



(B.15) 

(B.16) 
(B.17) 

(B.18) 



where we have used that nit = 175±10 GeV. Note that 5at is negligible and has the antiscreening 
sign (the screening of the t-quark loops in QED begins at q^ ^ m^, while in our case q^ = 
m\ < ni\). 

Finally, the W loop gives 



Ztx 



(3 + 4c2 



1 — v4(? — T arcsin — 
2c 



0.00050 



(B.19) 



The value of 5aw depends on gauge [^]; here we give the result of calculations in the 't Hooft- 
Feynman gauge. Traditionally, the definition of a takes into account the contributions of leptons 
and five light quarks only. The terms 5at and 5ay/ are taken into account in the electroweak 
radiative corrections. In our approach, these terms give the corrections 5iVi. In the same way 
the loops of not yet discovered heavy new particles ( "New Physics" ) should be accounted for. 
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Appendix C. 
How cxw{q^) and cxz{q^) 'crawl'. 

The effect of 'running' of electromagnetic coupling constants a{q'^) (logarithmic dependence 
of the effective charge on momentum transfer g^) is known for more than four decades) ||7^. In 
contrast to a(g^), the effective constants of W and Z bosons awi^'^) and az(g^) in the region 
< Q'^ ^ iTi'^z 'crawl' rather than run [Q. 

If we define the effective gauge coupling constant g'^{q'^) in terms of the bare charge g^ and 
the bare mass mo, and sum up the geometric series with the self-energy S(g^) inserted in the 
gauge boson propagator, this gives the expression 

Here m is the physical mass, and S(g^) contains the contribution of fermions only, since loops 
with W , Z and H bosons do not contain large logarithms in the region |g^| < m|. 

The bare coupling constant in the difference 1/ g'^{q'^) — 1/(7^(0) is eliminated, which gives a 
finite expression. The result is 

l/«z(9^) - l/ttz(0) = hzF{x) , where x = q^lm\ , (C.2) 

\law{(i) - l/ai^(0) = bwF{y) , where y = q^/m\ > (C-3) 

F(x) = — ^ log Ixl (C.4) 

1 — X 

If X 3> 1, equations ( p.2|) and ( |C.3|) define the logarithmic running of charges owing to 
leptons and quarks, and hz and hy^/ represent the contribution of fermions to the first coefficient 
of the Gell-Mann-Low function: 

hz = -^{iV„3[l + (l-^.Y]+iV.3[l + (-l + ^sy]+iVa2 + (l + (l-4.2f)]} , 
4o7r o o 

bw = T^ [GNg + 2Ni] , (C.5) 

where Nu^d,q,i are the numbers of quarks and leptons with masses that are considerably lower 
than y^. 

For g^ < m|, the numerical values of the coefficients bz^w are |]80| : 



6^-0.195 
bw ^ 0.239 



The massive propagator 2_ 2 in ( |C1| ) greatly suppresses the running of awiq ) and az{q 



Thus, according to ( |C.2|) and ( |C.3| ), the constant az(g^) grows by 0.85% from g^ = to 



2 2 

1/aziml) = 22.905 



l/az{ml) - 1/aziO) = -0.195 , (C.6) 
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and the constant awig"^) grows by 0.95%, 

l/aw{ml) = 28.74 

1/awiml) - l/awiO) = -0.272 , (C.7) 

while the electromagnetic constant a{q'^) increases by 6.34%: 

l/a(m|) - l/awiO) = 128.90 - 137.04 = -8.14 (C.8) 

With the accuracy indicated above, we can thus assume 

az{ml) ~ az{0) 
awim^z) — avi^(O). (C.9) 

At the same time, a{m\) differs greatly from q;(0); therefore the latter has no connection 
to the electroweak physics but only to the purely electromagnetic physics. 
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Appendix D. 

General expressions for one-loop corrections to 
hadronless observables. 

The bare quantities are marked by the subscript '0'. In the electroweak theory, three bare 
charges cq, /o and go that describe the interactions of 7, Z and W are related by a single 
constraint: 

(eoM)' + (W/o)' = 1 ■ (D.l) 

The bare masses of the vector bosons are defined by the bare vacuum expectation value of 
the higgs field 77: 

mzo = -foV , fTT-wo = -^9oV ■ (D.2) 

The fine structure constant a = e^/47r is related to the bare charge eo by the formula 



47r \ ^ c mi 



« = a(g2 = o) = ^(l-S;(0)-2-^^ , (D.3) 



where S'(0) = \im.q2^Q'E{q'^)/q'^. In the Feynman gauge S^^(O) ~ — (a/27r)(m^/cs)(l/e), 
where the dimension of spacetime is D = 4 — 2e. In the unitary gauge S^^(O) = 0. 

The simplest way to verify the presence of the term 2(s/c)S^2(0)/m| is by considering the 
interaction of a photon with the right-handed electron ept- Note that in this case there are no 
weak vertex corrections due to the PF-boson exchange. (Note also that the left-handed neutrino 
remains neutral even when loop corrections are taken into account, since the diagram with the 
7 — Z — uii^L interaction is compensated for by the vertex diagram with the W exchange). 

The relation between a = a(g^ = m|) and «o has the following form: 



a = ao 



ii^(m|) - s;(o) + s;(o) - 2^n^z(o) 



(D.4) 



where n^(g^) = S^(g^)/m|, n^^(g^) = S^^(g^)/m|, while E^ mean that contributions of 
W^-boson and t-quark are not accounted for. It is convenient to introduce in (p.4|) explicit 
expression for 6aw + So^t'- 

a = ao[l — n^(m|) — 2-11^^(0) — Saw — Sat] , (D.5) 



where in accordance with eq.( |B.3|) 



5a^ + 5a, = fl^(m|)-n^(m|) + S;(0)-S;(0) . (D.6) 

In the case of "New Physics" one should add to eq.( |D.5| ) the term Sa^p- Our first basic 



equation is eq.( p.5|) 



The second basic equation is: 

m| = m|o[l - Uziml)] = m\^Jcl[l - Iiz{ml)] . (D.7) 

A similar equation holds for m^: 

m^ =m^o[l -nty(m^)] , (D.8) 
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where Ui{q'^) = T.i{q^)/ml t = W,Z. 
Finally, the third basic equation is 



G 



9l 



^ A^m^ 



-[i + n^(o) + D] 



(D.9) 



wo 



where nvi/(0) = Svi/(0)/m^ conies from the propagator of W , while D is the contribution of 
the box and the vertex diagrams (minus the electromagnetic corrections to the four-fermion 
interaction) to the muon decay amplitude. According to Sirlin [l8l| , 



D 



a 



47rs2 ' 2s2 



where 



A^y = A(miy) 



■(6 + 
2 



7-4s2 



\ogc^ + A^w) 



+ log47r - 7 - log(m^//i2) 



(D.IO) 



(D.li: 



A-D 

Now we are able to express /o and qq in terms of a, G^, irtz and the loop corrections. From 
O, (PI ) and (^) we obtain: 



[nir 



fl = 4V2G^ml[l - IlwiO) + nziu^zJ 
From (PI) , dDj) and (pl^) we get: 



^] • 



(D.12) 



^0 

Co = "T 

Jo 



1 + 



2(c2 - S2) 



-2-n^z(0) - n^(m|) - 5aw - Sat + ^z{m\) - 0^/(0) - D) 



(D.13) 

The next step is to express my/ /mz, Qa and gy through c, s and loop corrections. Let us 
start with mw/'m-z- From ( |D.8| ) and ( |D.7] ) we get: 



"^wMz = Co[l - -Uwiml^) + -Uziml)] . 



(D.14) 



Substituting Cq given by ( |D.13| ) we obtain: 



mw 
mz 



c + 



cs 



2(c2 - s^) 



-[nz(m|) - Ilwim''^)] + Hwiml,) - 0^.(0)- 



n^(m 



^; — 2-11^2(0) — D — 6aw — Sai 



(D.15) 



In order to obtain expression for g^ we should recall that it is proportional to /o and take 
into account the Z boson wave function renormalization and Zll vertex loop correction: 



9A = -l- \[^z{ml) - nH.(0) -D- Y!^{m\) - ScsFa] 



(D.16) 



where Fa originates from the vertex correction. 

The last quantity is the ratio gv /qa- One-loop corrections come from Sq = 1 — Cg (eq.( p.l3"[ )), 
from vector and axial Zll vertices and from Z — > 7 transition which contributes to gv only: 



^ = 1 

9A 



As' 



A&s 



,2„2 



C^ — S^ " c 



[2-n^^(0) + n^(m|) + 5aw + 5at 



- Uziml) + Uw{0) +D]- AcsFv + AcsFa{1 - 4s') - 



— 4csn^z(m 



(D.17) 
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Formulas ( p. 151 ), ( p.l6| ) and ( |D.17| ) are derived in this Appendix according to the "Five 
Steps" procedure described in Section 4.6. They describe finite one-loop corrections to hadron- 
less observables. 

It is easy to evaluate the contribution of t-quark to physical observables in the approximation 
~ awf^t- In this approximation Uw^m'^) = ll]y{0), Ilz{jn%) = n^(0), 11^(0) — 11^^(0) = 
Sa/lGns'^c^ t and from eq.( p.l5| )-( |D.17| ) we get: 



mw 
mz 



c + 



Sac 



327r(c2 - s2)s2 



9A 



t) 



(D.18) 
(D.19) 

The corrections proportional to m^ were first pointed out by Veltman |3^, who emphasized 



1 3a; 

'r ^ 327rs2c2 



9A 



^ ~ 1 _ 4^2 + 



3a 



the appearance of such corrections for a large difference rn^ — ml which violates the isotopic 
symmetry. In this review the coefficients in front of the factors t in equations ( p.l8| )-( p.20| ) are 
used as coefficients for normalized radiative corrections Vi (see Sections 4.2 and 4.3). 



50 



Appendix E. 
Radiators Raq and Rvq- 

For decays to light quarks q = u,d, s, we neglect the quark masses and take into account the 
gluon exchanges in the final state up to terms ~ a^ [Q - |^, and also one-photon exchange 



in the final state and the interference of the photon and the gluon exchanges [gg]. These 
corrections are slightly different for the vector and the axial channels. 
For decays to quarks we have 

r, = r(Z -. qq) = 12[gl^RA, + gl^^Ry^n (E.l) 

where the factors Ra,v are responsible for the interaction in the final state. According to |p2| 



8l: 



Rv, = i + - + -g---gj-- + 

TT 4 ^ 7T 4 ^ IT IT 



s\2 



t IT 

- 12.77(^f + 12^^5.^ (E.2) 

IT rriz ^ 

Ra, = Rv, - {2T,,)[h{t)(^f + h{t){-f]- 

TT TT 



where rhq is the running quark mass (see below), 

<5„^ = l + 8.7(^)+45.15(^)^ (E.4) 

7r TT 

-^L = 1 + 3.67(^) + (11.29 - logt)(^)^ (E.5) 

TT TT 

2 8 logt 

^«- = 81 + ^' ^^-^^ 

/.(t) ^ -3.083 -logt+^ + ^, (E.7) 

/3(i) = -15.988- 3.722 logt+1.9171og^t , (E.8) 

t = ml /ml ■ 

Terms of the order of (d^/Tr)^ caused by the diagrams with three gluons in intermediate state 
were calculated in [^. For Ryq they are numerically very small, ~ 10~^; for this reason, we 
dropped them from formula ( [E.2| ). 
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For the Z -^ bb decay, the 6-quark mass is not neghgible; it reduces Tb by about 1 MeV 
(~ 0.5%). The gluon corrections result in a replacement of the pole mass ruh ~ 4.7 GeV by the 
running mass at q^ = m\ '■ ^^b -^ nii,{mz)- We express rhh{mz) in terms of rrib, as{fnz) and 



asirrib) using standard two- loop equations in the MS scheme (see p6|). 

For the Z ^ cc decay, the running mass rhcirriz) is of the order of 0.8 GeV and the 
corresponding contribution to r^ is of the order of 0.05 MeV. We have included this tiny term 
in the LEPTOP code, since it is taken into account in other codes (see, for example, @). 

We need to remark in connection with Fc that the term l2{t), given by equation ( [E.?] ), 
contains interference terms ~ (ds/vr)^. These terms are related to three types of final states: 
one quark pair, a quark pair and a gluon, two quark pairs. This last contribution comes to 
about 5% of I2 and is below the currently achievable experimental accuracy. Nevertheless, in 
principle these terms require special consideration, especially if these quark pairs are of different 
flavors, for example, bbcc. Such mixed quark pairs must be discussed separately. 

Note that d^ stands for the strong interaction constant in the MS subtraction scheme, with 
/i^ = m|. 
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Appendix F. 
a^t^ corrections from reducible diagrams. 



In [^ when deriving equations for physical observables we systematically took into account 
corrections which contained first power of polarization operators and neglected terms ~ {Ilw,z)'^- 
This procedure was correct at one loop, but since Ilw,z contain terms of the order of awt we 
evidently lost a^t^ terms. To restore them let us repeat the procedure implemented in |Q 
this time taking squares of IIi^ and Hz (reducible two- loop diagrams) into account. 

Our starting point are three basic equations for quantities rriz, G^ and a = a{rn\). Since 
terms ~ awt are absent in 11^, Il^z and D functions, we will not consider these functions in 
this Section. Equation for niz is the same as in Appendix D, eqs.( p.7[ ), ( p.2| ): 

^l = 7/oV[l-nz(m|)] , (F.l) 



4- 



while for G^ we have: 



and we keep nvi/(0) in the denominator to avoid loosing of the 11^(0) term (compare with 
eq.( |D.9| )). From these two equations we get: 

where we use equality Ilz{m\) = 11^(0) which is valid for the leading term ~ awt. 

Finally dividing the equation for the running electromagnetic coupling constant, which in 
our approximation is simply: 

e\ml) = Ana = glil - ^) , (F.4) 

Jo 

4(1 - 4) = ^^[1-^] , (F.5) 

/o Jo V2Gf,m^z 

i-n^(o) ^ nz{o)-Uw{o) 

- l-UwiO) l-UwiO) ■ ^ - ' 

Considering S as a small parameter and solving equation ( |1''.5| ) perturbatively, we get: 

I = c'[l + ^5 - -/^,6'] , (F.7) 

Jo c^ — s^ [c' — s^y 

where we keep terms linear and quadratic in 6. For definitions of c and s see eq.(^3D. 

Our next step should be the calculation of the 6'^ corrections to the functions Vi. But first 
let us discuss the expression for 6 as given by eq. (|F.6|) which contains factor 1 — nvK(O) in 
denominator. At one loop, corrections proportional to 6 appear in physical observables. They 
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by ( [b'.3[ ), we obtain: 



should be carefully calculated in order not to induce extra a^t^ terms. Fortunately this can 
be done straightforwardly, using the following chains of equalities: 



nz(0) = 7^/o'Pz(0)//o'] = 4V^G,(1 - UwiOmziO)//^] , 
JoV 



l^wiO) 



gW 



gl[i:w{Q)/9l] = 4%/2G,,(l - n^(0))[Ew^(0)^ 



0^] , 



(F.8) 



(F.9) 



where expressions in square brackets contain self-energies without coupling constants (S^/Zq 
and Yiw / Qq respectively) and eq.( |l:''.2| ) is used to express rj through G^. Substituting eq. 
and eq.( [F.9| ) into eq. (|F.6| ) we get: 



rrir. 



3a rrit .2 



(F.IO) 



Now everything is prepared for the calculation of 6"^ corrections to physical observables. Let 
us start from the H^-boson mass. For the ratio of the squares of vector boson masses we have: 



ml ~ fi l-Uziml) ■ 
Taking the ratio of bare coupling constants from equation (|F.7|) we get: 



:f.ii) 



mw 
mz 

It is easy to see that: 



T^wiO) 



[1 + 



\ 1 - UziO) '^ ' 2(c2 - s2 



-6 + 



^s^cH\ 



(C2 



.2)3 3 



(F.12) 



i-nH.(o) 
\ i-n^(o) 



i-nz(o) 
i-nw-(o) 



yr^i 



:F.13) 



Resulting formula for the correction to the ratio mw /rnz is presented in Section 7.1. 
The next step is the correction to axial coupling of Z boson to charged leptons. Axial 
coupling is proportional to /o, and from eq.( |l:''.3| ) and eq.( [b'.13| ) we immediately obtain: 



The final formula for the correction to qai is presented in Section 7.1. 

For the ratio of vector to axial constants in our approximation we have: 



(F.14) 



9vil9 



Al 



l-4s 



l-4fl- 



/o 



(F.15) 



The expression for the correction to gvi/OAi through physical parameters is presented in 
Section 7.1 as well. 
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Appendix G. 
Oblique corrections from new generations and SUSY. 

In this Appendix we collect analytical formulas for different oblique corrections. 
For the degenerate case the contribution of additional quark and lepton to 5^1^ are given 
by ( m): 

5'Vm = ^N,{[{l-l)F{l)-{l-l/c')F{l/c')]+2s'[{l-l/c')F{l/c')- 

- {l + 2l)F{l)]+As\Ql + Ql)[{l + 2l)F{l)-^-]^ (G.l) 



5'Va = -N,{[l-l + i6l'-3l)Fil)] + [As\Ql + Ql)- 

- 2s^][2l + l~l2l^F{l)]}/{l-Al) , (G.2) 

5'Vn = -^iV, |3/F(/) - As'c'iQl + Ql) [(1 + 2/)F(/) -\]\ , (G.3) 

where N^ = 3, Qjj = 2/3, Qd = ~l/3 for quark doublet; Nc = 1, Qu = 0, Qd = —1 for lepton 
doublet; 

i = m^Q/m\ fo^ quarks , £ = m\/m\ for leptons , 

and the function F{1) is defined in Appendix B, eqs.(B8), (B.IO). 

For different up- and down- quark (and lepton) masses analytical expressions for 5^Vi are 
given in by 

1 64 Ifi 2 

-5'V^ = C-s'-^s')[{l + 2u)F{u) + {l + 2d)F{d)--] + 

+ ^[(1 - u)F{u) + (1 - d)F{d) - ^] + ^^[« + ^ - f^ log ^] + 
8, s"^ .,u — d ^ u , ,, ,, u + d.u + d , u 4, 

- (2c'-i<-<i-^^^^^)FK,,mJ,m|,)l ; (G.4) 



4 r 16 
3^Va = -\{—s^-As^-l)[2uF{u)-{l + 2u)F'{u)+2dF{d)-{l + 2d)F'{d)] + 
9 13 



2vd 7/ 1 

+ 3[u + d log -F'{u)-F\d)]\ ■ (G.5) 

u — d d J 

U'Vn = -l[uF{u)+dF{d) + -^^\ogl-^^] + ^^s'c\^^ , 

(G.6) 
where n is the number of generations and rriN = f^u, rriE = rno, u = ra^/m^, d = m\)/m\\ 
F'{u) = —u{d/du)F{u) and F(s,m^,m2) is defined in Appendix B, eq.(B.6). 
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The formulae that describe the enhanced SUSY corrections to the functions Vj have the 



following form ||72|| : 



^susyVa = -^[clg{mi,mi) + s^c/(m2,m-J - clslg{mi,m2)] , (G.7) 

1 777 777 1 

SI§syVr = S^s§syVa + ^Y,[cllog{-^) + sl\og{J)] - -clslhim„m,) , (G.8) 

O III- 1 1 It O 



5's^syV^ = Sl§sYVA + lY,s^[cl\ogi^) + sl\ogi^,)] + 

O TTh~ TTh~ 

o 

+ ^-^^[c^/i(mi,m^) + s2/i(m2,m^)]-^/i(mi,m2) , (G.9) 

where 

2 2 2 

^(mi,m2) =m? + m2-2^Ji^log(^) , (G.IO) 

nil ~" ""^2 ^^2 

and Yl = Qt + <5b = 1/3 is the hypercharge of the left doublet. 
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Appendix H. 
Other parametrizations of radiative corrections. 

Here we will present formulas which connect our functions Vi with two other sets of param- 
eters widely used in the literature to describe electroweak radiative corrections. All formulas 
of this Appendix are valid at one electroweak loop approximation. 

A set of three parameters ei, 62, £3 has been suggested by Altarelli, Barbieri and Jadach 



87|| for phenomenological analysis of New Physics: 

ei = Ap , (H.l 



s^ 



e2 = c^Ap+- -Arw~2s^Ak' , (H.2 

c^ — s^ 

£3 = c^Ap + (c^ - s^) AA;' , (H.3 
where Ap describes the correction to qa, Ak' to gy and Ar^y to mw /mz'- 

gA = -\{l + \Ap) , (H.4 

gv/gA = l-^s'{l + Ak') , (H.5 

mw/mz = c[l-s^Arw/2{c^-s^)] . (H.6 
By comparing these definitions with the definitions of Va, Vr and Vm we obtain: 

. 3a Va ^^^ 

^k' = -i^ ,^, '" ^, , (H.8 



Hence: 



^H 


IQtt s'-' 


1^2 ' 


{ 


3a 


Vr 




167r (c2 


- S2)s2 


Arw 


3a Vm 
167r s^ 


Bi = 


3a 

1 /^ 9-3 


-Va , 



(H.9 



(H.IO 

±U/I S'C- 

e2 = TTTT^^^T-^liVA - Vm) - 2s'iVA - Vr)] , (H.ll 

^3 = -^,{Va - Vr) . (H.12 

As is evident from the last two formulas, the virtue of £2 and £3 is that they do not contain the 
term t. So, at the time when top-quark mass was not measured at Tevatron the corresponding 
uncertainties in 62 and £3 were diminished. 

Another set of parameters, S, T, U was introduced a few years earlier by Peskin and 
Takeuchi [BH] . These parameters were proposed to describe only the so-called oblique corrections 



due to the physics beyond the Standard Model. Using the definitions of S, T, U from |88] and 
designating New Physics contributions to e^ as fej we obtain: 

fei = aT , (H.13) 
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fe2 = -aU/As^ , (H.14) 

6e3 = aS/4s'^ . (H.15) 



From (plOD - ( W^ ) we get: 

3 



47r(c"^ — s^j 

S = ^{5Va-5Vr) , (H.18) 

47r 



where 5Vi are New Physics contributions to Vi. 



According to |g8 



S = 167r[S;,3(0) - S'3q(0)] = 167r[S'^(0) - S^O)] , (H.19) 

T=-^[Sn(0)-S33(0)] , (H.20) 

f/ = 167r[S;i(0)-S'3(0)] , (H.21) 

where S'(0) = dT,{q'^) /dq^ \g2^Q and S's are defined by the corresponding currents (isotopic, 1 
and 3, and electromagnetic, Q, vector, V, and axial. A) with coupling constants being extracted. 
Thus S characterizes the degree of chiral symmetry breaking, while both T and U that of 
isotopic symmetry. Note that in eqs. ([H.19|) - ( |H.21| ) only the contribution of New Physics 



should be considered. Since new particles should be heavy it is reasonable to take into account 
only values of self-energies at g^ = and their first derivatives (higher derivatives are power 
suppressed) |5^. Altogether we have eight parameters (T,^r^r(^0), T,zz{0), T,^z{0), ^^7(0), 
^Wvk(O)' ^zz(0)5 ^^7(0)5 ^77(0)) 5 two of which are equal zero (S^^(O) and T,^z{0)), while 
three combinations can be absorbed in the definition of a, G^ and mz- The remaining three 
combinations enter S, T and U (or fei, z = 1, 2, 3). 
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Figure Captions 

Figure 1: Muon decay in the tree approximation. 

Figure 2: The Z boson as a resonance in e^e~ annihilations. 

Figure 3: Photon polarization of the vacuum, resulting in the logarithmic running of the 



2 



electromagnetic charge e and the 'fine structure constant' a = 1^, as a function of g , where q 
is the 4-momentum of the photon (a). Some of the diagrams that contribute to the self-energy 
of the W^-boson (b)-(g). Some of the diagrams that contribute to the self-energy of the Z-boson 
(h)-(n). Some of the diagrams that contribute to the Z <-> 7 transition (o)-(r). 

Figure 4: Vertex triangular diagrams in the Z ^ II decay (a), (b), (c). Loops that renor- 
malize the lepton wavefunctions in the Z ^ II decay. (Of course, antilepton have similar loops.) 
(d), (e). Types of diagrams that renormalize the Z boson wavefunction in the Z ^ II decay 
(f), (g). Virtual particles in the loops are those that we were discussing above. 

Figure 5: Virtual t-quarks (a) and W bosons (b), (c) in the photon polarization of the 
vacuum. 

Figure 6: Gluon corrections to the electroweak quark loop of the Z-boson self-energy. 

Figure 7. Kn as a function of mj for three values of mu- The dotted parabola corresponds 
to Veltman approximation: Vm = t. Solid horizontal line traces the experimental value of Vm 
while the dashed horizontal lines give its upper and lower limits at the la level. 

Figure 8. Va as a function of mt- The dotted parabola corresponds to Veltman approx- 
imation: Va = t. Solid horizontal line traces the experimental value of Va while the dashed 
horizontal lines give its upper and lower limits at the la level. 

Figure 9: Vr as a function of rrit. The dotted parabola corresponds to Veltman approxi- 
mation: Vr = t. Solid horizontal line traces the experimental value of Vr while the dashed 
horizontal lines give its upper and lower limits at the la level. 

Figure 10: rrit — mu exclusion plots with assumptions of: a) irit = 150(5) GeV; b) rrit = 
175(5) GeV; c) rrit = 200 ± 5 GeV. 

Figure 11: The vertex electroweak diagrams involving t-quark and contributing to the 
Z ^ bb decay. Diagram (b) represents gluon corrections to diagram (a). 

Figure 12: Some of Feynman diagrams that give a^rt'^ corrections. 

Figure 13: x^ ''^s- ^^h curves. 

Figure 14: The 2-dimensional exclusion plot for the case of A^ extra generations and for the 



choice rriD = 130 GeV - the lowest allowed value for new quark mass from Tevatron search [27 



using niH > 90 GeV at 95 % C.L. from LEP-2 |p9[. Little cross corresponds to x^ minimum 



lines show one sigma, two sigma, etc allowed domains. 

Figure 15: Contribution of t and b squarks into W- and Z bosons self-energy. 
Figure 16: Values of SVa, SVr and 6Vm at m^ = 200 GeV. 
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Figure 1: Muon decay in the tree approximation 




Figure 2: The Z-boson as a resonance in e^e annihilations. 
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Figure 3: Photon polarization of the vacuum, resulting in the logarithmic running of the 
electromagnetic charge e and the 'fine structure constant' a = 1^, as a function of g^, where q 
is the 4-momentum of the photon (a). Some of the diagrams that contribute to the self-energy 
of the PF-boson (b)-(g). Some of the diagrams that contribute to the self-energy of the Z-boson 
(h)-(n). Some of the diagrams that contribute to the Z ^^ 7 transition (o)-(r). 
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Figure 4: Vertex triangular diagrams in the Z ^ II decay (a), (b), (c). Loops that renormahze 
the lepton wavefunctions in the Z ^ II decay. (Of course, antilepton have similar loops.) (d), 
(e). Types of diagrams that renormalize the Z boson wavefunction in the Z ^ II decay (f), 
(g). Virtual particles in the loops are those that we were discussing above. 
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Figure 5: Virtual t-quarks (a) and W bosons (b), (c) in the photon polarization of the vacuum. 



69 



-^^{jy^ ^ 



a) 




Figure 6: Gluon corrections to the electroweak quark loop of the Z-boson self-energy. 



70 



> 




00 125 150 175 200 225 250 

mtop (GeV) 



Figure 7: Vm as a function of rrit for three values of itlh- The dotted parabola corresponds 
to Veltman approximation: Vm = t. Solid horizontal line traces the experimental value of Vm 
while the dashed horizontal lines give its upper and lower limits at the la level. 
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Figure 8: V^ as a function of rrit. The dotted parabola corresponds to Veltman approximation: 
Va = t. Solid horizontal line traces the experimental value of Va while the dashed horizontal 
lines give its upper and lower limits at the la level. 



72 



> 




250 



Figure 9: Vr as a function of mt. The dotted parabola corresponds to Veltman approximation: 
Vr = t. Solid horizontal line traces the experimental value of Vr while the dashed horizontal 
lines give its upper and lower limits at the la level. 
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Figure 10: rrit — mu exclusion plots with assumptions of: a)mt = 150(5) GeV; b) m^ = 175(5) 
GeV; c)mt = 200(5) GeV. 



74 





a) 



b) 



Figure 11: The vertex electroweak diagrams involving t-quark and contributing to the Z 
decay. Diagram (b) represents gluon corrections to diagram (a). 
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Figure 12: Some of Feynman diagrams that give a^t corrections 
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Figure 13: x^ vs. tuh curves. 
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Figure 14: The 2- dimensional exclusion plot for the case of N extra generations and for the 
choice mo = 130 GeV - the lowest allowed value for new quark mass from Tevatron search (p7[|), 
using rriH > 90 GeV at 95 % C.L. from LEP-2 |^. Little cross corresponds to x^ minimum; 
lines show one sigma, two sigma, etc allowed domains. 




Figure 15: Contribution of t and b squarks into W- and Z bosons self-energy. 
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Figure 16: Values of the 6Va, SVji and 6Vm at mr = 200 GeV. 
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